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DYNAMIC SUBSTRUCTURING BY THE BOUNDARY FLEXIBILITY 
VECTOR METHOD OF COMPONENT MODE SYNTHESIS 


Abstract 

by 

AYMAN AHMED ABDALLAH 


Component mode synthesis (CMS) is a method of dynamic 
analysis, for structures having a large number of degrees of 
freedom (D.O.F.). These structures often required lengthy 
computer CPU time and large computer memory resources, if solved 
directly by the finite-element method (FEM). In CMS, the 
structure is divided into independent components in which the 
D.O.F. are defined by a set of generalized coordinates defined by 
displacement shapes. The number of the generalized coordinates 
are much less than the original number of physical D.O.F. , in the 
component. The displacement shapes are used to transform the 
component property matrices and any applied external loads, to a 
reduced system of coordinates. Reduced system property matrices 
are assembled, and any type of dynamic analysis is carried out in 
the reduced coordinate system. Any obtained results are back 
transformed to the original component coordinate systems. In all 
conventional methods of CMS, the mode shapes used for components. 
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are dynamic mode shapes, supplemented by static deflected shapes. 
Historically, all the dynamic mode shapes used in conventional CMS 
are the natural modes (eigenvectors) of components. 

This work presents a new method of CMS, namely the boundary 
flexibility vector method of CMS. The method provides for the 
incorporation of a set of static Ritz vectors, referred to as 
boundary flexibility vectors, as a replacement and/or supplement 
to conventional eigenvectors, as displacement shapes for 
components. The generation of these vectors does not require the 
solution of a costly eigenvalue problem, as in the case of natural 
modes in conventional CMS, and hence a substantial saving in CPU 
time can be achieved. The boundary flexibility vectors are 
generated from flexibility (or stiffness) properties of 
components. The formulation presented is for both free and 
fixed-interface components, and for both the free and forced 
vibration problems. Free and forced vibration numerical examples 
are presented to verify the accuracy of the method and the saving 
in CPU time. Compared to conventional methods of CMS, the results 
indicate that by using the new method, more accurate results can 
be obtained with a substantial saving in CPU time. 


iii 



To my Mother, my Wife 



ACKNOWLEDGEMENTS 


I would like to express my deepest and sincere gratitude to my 
advisor Dr. Huckel bridge, for his guidance, advice, suggestions 
and friendship. Working with him was an enlightening experience 
for me. I would like to express my appreciation to Dr. D. 
Gasparini, for being an excellent teacher in all courses I have 
taken with him, thus providing me with an excellent foundation in 
my area of research and in other areas of great interest to me. 
Also I would like to thank him for being a member of my graduate 
committee. Special thanks are due to Dr. C. Lawrence, at 
NASA-Lewis Research Center, for his thoughtful suggestions and 
being a member of my graduate committee. I would like to 
acknowledge the members of the graduate committee, Dr. R. Mullen 
and Dr. R. Quinn for their valuable suggestions. I wish to 
express my gratitude to my wife, Nevine, for her understanding, 
love, support and patience in typing the text. 

The support of this research work by the NASA Lewis Research 
Center at Cleveland, Ohio, under grant NAG 3-707, is gratefully 
acknowledged. 


v 



Table of Contents 


page 

Abstract i i 

Dedication iv 

Acknowledgements v 

Table of Contents vi 

List of Figures ix 

List of Tables xiii 

Chapter 1. Introduction 1 

Chapter 2. Literature Review 5 

2.1. Introduction 5 

2.2. Review of Component Mode Synthesis :.... 5 

2.3. Review of Previous Work with Load-Dependent 

Ritz Vectors 10 

Chapter 3. Boundary Flexibility Vector Method of 

Component Mode Synthesis 19 

3. 1. Introduction 19 

3.2. Fixed-Interface Boundary Flexibility Vector 

Method of CMS 20 

3.2.1. Boundary-Flexibility Vector 

Formulation 20 

3.2.2. Hybrid Boundary Flexibility Vector/ 

Conventional CMS Formulation 25 

3.3. Free-Interface Boundary Flexibility Vector 

Method of CMS 27 

3.3.1. Boundary Flexibility Vector 

vi 



Formulation 28 

3. 3. 1.1. Components Having No Rigid-Body 

Modes 28 

3.3. 1.2. Components Having Rigid-Body 

Modes 38 

3.3.2. Hybrid Boundary Flexibility Vector/ 

Conventional CMS Formulation 41 

3.4. Computational Effort for Boundary Flexibility 

Vectors vs Normal Modes of Vibration 43 

3.5. Summary 46 

Chapter 4. Results of Numerical Problems in Free 

Vibration : . . . . 48 

4. 1. Introduction 48 

4.2. Sample Problem One 48 

4.2.1. Solution by Fixed- Interface Methods 

of CMS 49 

4.2.2. Solution by Free-Interface Methods 

of CMS 51 

4.3. Sample Problem Two 53 

4.4. Sample Problem Three 55 

Chapter S. Forced Vibration Response by The Boundary 

Flexibility Vector Method of CMS 74 

5.1. Introduction 74 

5.2. The Application of The Boundary Flexibility 
Method of CMS in Forced Vibration 

Problems 74 

vii 



5.3. A Comment on Wilson’s Load-Dependent Ritz 

Vectors 82 

5.4. Sample Problem 86 

Chapter 6. Summary and Conclusion 115 

References 118 

Appendix A Comparison of The Number of Operations For 

Generating The Boundary Flexibility Vectors 

Versus Eigenvectors 120 

Appendix B Flexibility Matrix of Components In Terms of 

Stiffness Submatrices 124 

Appendix C Proportional Damping Matrix 126 


viii 



List of Figures 


Figure p age 

2.1 A Structural Component 18 

4. 1 Sample Problem One 58 

4.2 Component Constraint Modes 59 

4.3 Component Attachment Modes 60 

4.4 Variation of Error in o> with K for Fixed 

3 

Methods 61 

4.5 Variation of Error in w with K for Fixed 

4 

Methods 62 

4.6 Variation of Error in o) with K for Fixed 

s 

Methods 63 

4.7 Variation of Error in o) with K for Free 

3 

Methods 64 

4.8 Variation of Error in u with K for Free 

4 

Methods 65 

4.9 Variation of Error in w with K for Free 

5 

Methods 66 

4. 10 Sample Problem Two 67 

4.11 Sample Problem Three 68 

4. 12 Component Number One for Sample Problem 

Three 69 

4. 13 % Saving in CPU Time for Generating Equal 

Numbers of Boundary Flexibility Vectors vs. 

Eigenvectors 70 

ix 



4.14 % Saving in CPU Time for Generating The Boundary 

Flexibility Vectors vs. One Eigenvector 71 

5. 1 Load for Sample Problem 90 

5.2 Response at Node 13 in The Z-Direction for The Case 
of Load Acting at Node 31 in The Z-Direction 

(Low Damping) 91 

5.3 Response at Node 31 in The Z-Direction for The Case 
of Load Acting at Node 31 in The Z-Direction 

(Low Damping) 92 

5.4 Response at Node 112 in The Z-Direction for The Case 
of Load Acting at Node 31 in The Z-Direction 

(Low Damping) 93 

5.5 Response at Node 13 in The Z-Direction for The Case 
of Load Acting at Node 31 in The Z-Direction 

(High Damping) 94 

5.6 Response at Node 31 in The Z-Direction for The Case 
of Load Acting at Node 31 in The Z-Direction 

(High Damping) 95 

5.7 Response at Node 112 in The Z-Direction for The Case 
of Load Acting at Node 31 in The Z-Direction 

(High. Damping) 96 

5.8 Response at Node 10 in The Z-Direction for The Case 
of Load Acting at Node 10 in The Z-Direction 

(Low Damping) 97 

5.9 Response at Node 84 in The Z-Direction for The Case 
of Load Acting at Node 10 in The Z-Direction 


x 



(Low Damping) 98 

5.10 Response at Node 112 in The 2-Direction for The Case 
of Load Acting at Node 10 in The Z-Direction 

(Low Damping) gg 

5.11 Response at Node 10 in The Z-Direction for The Case 
of Load Acting at Node 10 in The Z-Direction 

(High Damping) 100 

5. 12 Response at Node 84 in The Z-Direction for The Case 

of Load Acting at Node 10 in The Z-Direction 

(High Damping) 101 

5.13 Response at Node 112 in The Z-Direction for The Case 

of Load Acting at Node 10 in The Z-Direction 

(High Damping) 102 

5. 14 Response at Node 16 in The Y-Direction for The Case 

of Load Acting at Node 16 in The Y-Direction 

(Low Damping) 103 

5. 15 Response at Node 87 in The Y-Direction for The Case 

of Load Acting at Node 16 in The Y-Direction 
(Low Damping) 104 

5.16 Response at Node 99 in The Y-Direction for The Case 
of Load Acting at Node 16 in The Y-Direction 

(Low Damping) 105 

5.17 Response at Node 16 in The Y-Direction for The Case 
of Load Acting at Node 16 in The Y-Direction 

(High Damping) 106 

5. 18 Response at Node 87 in The Y-Direction for The Case 


xi 



107 


5. 19 


of Load Acting at Node 16 in The Y-Direction 

(High Damping) 

Response at Node 99 in The Y-Direction for The Case 
of Load Acting at Node 16 in The Y-Direction 
(High Damping) 


108 


xii 



Chapter 1 


Introduction 

Since the late fifties and early sixties, there has been 
rapid advancement in the field of structural dynamics. The 
primary cause of this advancement was the introduction of digital 
computers and the associated analytical software which made the 
application of numerical methods to dynamic analysis problems more 
feasible. The most powerful of these numerical techniques was the 
finite element method, which c an be considered the basis for a new 
era in the field of mathematics (more precisely in the solution of 
partial differential equations), and consequently many fields of 
science, including structural dynamics. 

Previous to the finite element method, there existed other 
approximate methods of analysis, such as Rayleigh's method in 
dynamic analysis. The main disadvantage of this type of analysis 
is that it requires the assumption of reasonably accurate mode 
shapes for the structure, which is not always readily 
accomplished, especially in large and complex structures. 
Compared to these previous methods of analysis, the finite element 
method is more accurate, more direct (no mode shapes need to be 
assumed) and can be applied to a structure of any degree of 
complexity. 

Following the finite element method, a new modeling 
technique, namely component mode synthesis (CMS), was introduced 
by Hurty [1] in the mid-sixties. Component mode synthesis is a 
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method of dynamic substructuring, in which the structure is 
divided into independent components or substructures. The main 
motivation for introducing the method arose in the dynamic 
analysis of large structures, having a large number of physical 
degrees of freedom. These large structures usually required 
lengthy computer CPU time and large computer memory resources, 
when solved by finite elements directly. By using CMS, savings in 
computer time and memory can be attained. Moreover, other 
advantages can also be attained through the use of dynamic 
substructuring such as: 

1- Independent design and analysis efforts for various components 
of a structure. 

2- Reducing effort if system can be divided into identical 
repetitive components having the same constraints and 
displacement modes. 

3- The incorporation of experimental and analytical data can be 
achieved, for the characterization of components. 

4- Reducing computational cost of a reanalysis, in case not all 
components are modified. 

The saving in CPU time and computer memory is achieved by 
reducing the number of degrees of freedom, associated with each 
component. This is accomplished by defining the displacements of 
each component by a set of generalized coordinates, consisting of 
the amplitudes of a corresponding set of mode shapes. In all 
conventional CMS methods, the mode shapes used are dynamic mode 
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shapes, supplemented by static deflected shapes. All the dynamic 
modes historically used in CMS, are natural inodes of vibration 
(eigenvectors) for the component. 

This work presents a new method of CMS, namely the boundary 
flexibility vector method of CMS, which is based upon utilizing a 
set of static Ritz vectors, as displacement shapes for components. 
The generation of these vectors, referred to herein as “boundary 
flexibility vectors", does not require the solution of a costly 
eigenvalue problem, as in the case of natural modes in traditional 
CMS. Thus, a substantial saving in computer time is obtained. 
The presented method is general and can be utilized for any type 
of dynamic analysis. 

This work is divided into six chapters. Chapter one is an 
introduction to the work. Chapter two is a brief review of 
traditional methods of CMS, and previous work in the application 
of static Ritz vectors in structural dynamics. The formulation of 
the boundary flexibility method of CMS is presented in chapter 
three. The presented formulation is for both fixed and 
free- interface components. Also, a hybrid boundary flexibility 
and traditional CMS formulation is presented. A comparison of the 
number of operations required to generate the boundary flexibility 
vectors versus eigenvectors is given at the end of the chapter. 
In chapter four, small and large free vibration problems are 
utilized to verify the accuracy and CPU time saving of the new 
method compared to finite element representations and traditional 
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methods of CMS. Chapter five presents the formulation of the 
boundary flexibility method of CMS for the forced vibration 
problem. An example is solved numerically by several methods and 
for different loading conditions. The transient response is 
compared for all methods. Chapter six includes a work summary, 
conclusions and ideas for future work in this field. 



Chapter 2 


Literature Review 

2 . 1) Introduction : 

This chapter is divided into two main and independent parts. 
The first part is a very brief review of component mode synthesis. 
The main steps of CMS, classif icat ion of components and type of 
displacement shapes are introduced in this part. Part two is a 
detailed review of previous work, in the application of load 
dependent Ritz vectors in structural dynamics. 

2 .2) Review of Component Mode Synthesis : 

As mentioned before in chapter one, component mode synthesis 
(CMS) was introduced by Hurty [1] in 1965. Since then, several 
different variations of the method were introduced. However, the 
main idea of reducing the number of degrees of freedom used to 
characterize components, through use of a truncated set of mode 
shapes and generalized coordinates, is common to all CMS methods. 

A detailed review of CMS methods was given in a previous work [2] 
by the author. The main steps of conventional CMS, as given in 
reference [2], are summarized as follows: 

1- Divide the system into independently characterized components. 
Each component has a vector {u} of total physical degrees of 
freedom (D.O.F. ). (See Figure (1)). A component can be connected 
to other components or to the ground through a support, as shown 
in Figure (1). The common boundary between components is the 
interface. Let {u } be the vector of D.O.F. contained in the 
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interfaces of the component. Thus, {u} is partitioned into 


where 


<V 


<u> 



^ i ' 

= vector of internal D.O.F. 

= complement of {u } in {u>. 


in component. 


( 2 . 1 ) 


2- Define the displacement of each component by a set of mode 
shapes and generalized coordinates. The generalized coordinates 
can be taken as the amplitudes of these mode shapes. The mode 
shapes used contain static displacement shapes, as well as a 
truncated set of natural modes of vibration. Usually a very small 
number of static and dynamic displacement shapes is required, to 
obtain a good representation of the component in the new 
coordinate system. Thus, the number of the new set of generalized 
coordinates is generally much smaller than the original number of 
physical D.O.F. , in the component. 

3- For each component, transform the property matrices (mass, 
stiffness and damping), from the physical coordinate system to the 
new reduced system of generalized coordinates. 

4- Couple all the components by enforcing displacement 
compat ibi 1 ity requirements at the interfaces between components, 
thereby assembling reduced system property matrices. 

5- Perform any type of dynamic analysis (eigenproblem solution, 
transient or steady-state response, etc. .. ) in the reduced 
coordinate system of generalized coordinates. 

6- Any obtained results are back transformed to the original 
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system of physical D.O.F. 

For all variations of C.M. S. , within the individual 
characterization of components, interface D.O.F. car be considered 
either free or totally fixed. Accordingly, components can be 
classified into either fixed or free- interface components. For a 
detailed discussion and formulation of fixed- interface components 
see Hurty [1], Abdallah [2], Craig and Bampton [3] and Craig [4]. 
While for free-interface components see Abdallah [2], MacNeal [5], 
Rubin [6] and Martinez and Gregory [7], The first investigator to 
introduce fixed- interface components, was Hurty [1] in 1965, while 
free-interface components were introduced by Goldman [8] in 1969. 
The method of coupling components, introduced by Goldman [8], was 
rather complicated. Several alternative methods of coupling 

free-interface components were introduced by Hou [9], Dowell [10] 
and Abdallah [2]. In 1971, following Goldman's work, MacNeal [5] 
gave a better understanding and representation of free-interface 
components through the use of the idea of residual flexibility. 
In 1975, Rubin [6] expanded MacNeal 's work to free-interface 
components, having rigid body modes. 

It should be noted that the decision to characterize 
interface D.O.F. into free or fixed often depends on the type of 
displacement shapes available to the analyst, i.e., the fixed or 
free-interface type. It also depends on the judgment of the 

analyst as to what level of constraint each component imposes on 
the other. The displacement shapes used for the characterization 
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of components, may be obtained either experimentally or 
analytically. There are several different types of displacement 
shapes used in the literature of C. M. S. They are summarized in 
the following: For fixed-interface components: 

1- Fixed- interface normal modes of vibration [<& ]: They are the 

n 

natural modes of vibration (or the eigenvectors), obtained from 
solving the eigenvalue problem associated with the equation of 
motion of the component. They are computed assuming all the 
interface degrees of freedom fixed. Only a truncated set of 
normal modes [$^] (kept modes) are used. Usually the low 

frequency mode shapes are the ones used or kept. 

2- Constraint modes [$ J: A constraint mode {<t> } is defined as, 

c c 

the static displaced shape obtained by applying a unit 
displacement at one interface D.O.F. while totally constraining 
all other interface D.O.F. in the component. The number of 
constraint modes used for each component is equal to the number of 
interface D.O.F. in that component. The generalized coordinates 
associated with these constraint modes, are usually taken as the 
interface D.O.F. in components. Rigid-body modes car be included 
in these constraint modes. 

For free-interface components, the following displacement shapes 
sire used: 

1- Free-interface normal modes of vibration [$ ]: They are the 

n 

same as in fixed-interface components except they are computed 
assuming all interface D.O.F. free. 



9 


2- Rigid-body modes [$^]: They ape used for free- interface 

components having no constraints preventing them from rigid-body 
motion. 

3- Attachment modes [<f ]: They are defined on a subset {u > of the 

a a 

total physical degrees of freedom {u>. An attachment mode { <p > is 


defined as the static displaced shape obtained by imposing a unit 
force (or moment) on one D.O.F. of {u > with zero forces on all 


remaining D.O.F. in {u^}. According to Meirovitch [11], Bamford 
was the first investigator to introduce attachment modes in 1967. 
They were used to account for the effect of concentrated loads on 
unconstrained D.O.F. 

4- Residual flexibility modes ['£ ]: They are static displacement 

d 

shapes introduced by MacNeal [5] in 1971, to obtain a better 
representation for free-interface components. (See (5)). They 
were introduced to replace attachment modes as static mode shapes 
for the component. 

5- Inertia relief attachment modes [$ ]: They ape the static 
displacement shapes, for components having rigid-body modes, 
analogous to attachment modes ] in case of components having no 


rigid-body modes. They were introduced by Rubin [6] in 1975. 

6- Residual inertia relief attachment modes [$ ]: They are static 

r 

displacement shapes introduced by Rubin [6] in 1975. They were 
introduced to be used for free- interface components having 
rigid-body modes analogous to residual flexibility modes [# ], 

d 

used for components having no rigid-body modes. 
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2 .3) Review of Previous Work with Load-Dependent Ritz Vectors : 

The application of load-dependent Ritz vectors was introduced 
by Wilson et al. [12] in 1982. They were the first investigators 
to point out the usefulness of load-dependent Ritz vectors in 
dynamic analysis. The vectors were used to calculate transient 
dynamic response by Ritz vector superposition. The load-dependent 
Ritz vectors were used, instead of natural modes of vibration 
(eigenvectors), in a fashion equivalent to mode shapes. In all 
presented examples, it was shown that the superposition of 
load-dependent Ritz vectors yielded more accurate forced dynamic 
response results, with fewer number of vectors, than if the 
natural modes (eigenvectors) were used. Furthermore, the 
load-dependent Ritz vectors are generated with less computational 
effort than eigenvectors. The method used by Wilson et. al [12] 
is summarized in the following: 

Consider a damped system, having n degrees of freedom. Its 
equation of motion is: 

CM] iu} + [C] {u> + [K] {u> = {f(s)> g(t) (2.2) 

where [M] = mass matrix of system 

[C] = damping matrix of system 
[K] = stiffness matrix of system 
{u} = displacement vector 
{u> = velocity vector 
{u} = acceleration vector 

{f(s)} = vector of spatial distribution of 




11 


external load 

g(t) = time-dependent amplitude of {f(s)> 

The load-dependent Ritz vectors are generated by Wilson’s 
algorithm given in reference [12]. The main steps ape as follows: 

1- Given [M], [K] and {f(s)>. 

2- Solve for first vector {x^ and normalize with respect to the 
mass matrix to obtain {x^. 

[K] { x*} = {f(s)> 

{x : > = {x t > / ( {x^ 1 [M] {x : > ) 

* X * 

where {x^} is the transpose of {x^} 

3- Solve for additional vectors (i = 2,3,4 L) 

[K] {x*> = [M] {x i > x 

- Orthogonalize with respect to [M], with all previous ( i— 1 ) 
vectors 

•• • * p* 

{x t > = {x > - y ( c <x }) 

1 1 j=l J J 

where c^ = {x^} 7 [M] {x*> 

- Normalize vector {x i } with respect to [M] obtaining {x^ 

1/2 

{x t > = {x”} / ( {x*V [M] {x’*> ) 

4- Orthogonalize all obtained vectors with respect to the 

stiffness matrix [K] (optional step). This step is done as 


follows: 
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- Form [X] = matrix containing all previously calculated 

load-dependent Ritz vectors in its columns. Its 
dimension is (n x L) 

[A] = [X] T [K] [X] 

[I] = [X] T [M] [X] = identity matrix of order L 

- Solve the following reduced-order (L x L) eigenvalue problem. 

(Solve for (y^, i = 1 L) 

[A] ~ Wj [II J {y t > = {0} 

where {y^> = i th eigenvector 

2 . th . 

w i eigenvalue 

- Compute final orthogonal vectors [X°] from 

tx°] = [x] m 

o 

where [X ] = an (n x L) matrix containing final 
load-dependent Ritz vectors in its columns. 
All vectors are orthogonal with respect to 
the mass and stiffness matrices. 

[Y] » an (L x L) matrix containing the 

eigenvectors {y^} in its columns. 

According to Wilson's algorithm, load-dependent Ritz vectors 
can be defined as static displacement shapes generated from the 
spatial distribution of the external loading applied on the 

system. It can be seen that the initial load-dependent Ritz 
vector is the static displacement of the system, due to the 
applied load, while all subsequent vectors are obtained from the 
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inertial loading of the response. Thus, it is ensured that all 
the generated load-dependent Ritz vectors, are excited by the 
applied external load pattern, and hence, will contribute to the 
response. While in the case where natural modes of vibration 
(eigenvectors) are used, there could be mode shapes which 
contribute little or nothing to the response (as the case of a 
symmetric spatial distribution of load and an anti-symmetric mode 
shape) even if the natural frequency of the mode shape is well 
represented in the loading frequency content. 

After obtaining the load-dependent Ritz vectors [X°], the 
displacement vector (u) is approximated by 

L 

(u) as V {X°> 2 (t) (2.3) 

1 = 1 1 1 

where 2^(t) = t ime-dependent generalized coordinates 

associated with load-dependent vector i. 

The mass, stiffness and damping matrices are transformed to a 
reduced set of coordinates by: 

CM*] = [X°] T CM] [X] 

= identity matrix tl] of order L 

[K*] = [X°] T [K] [X] 

[C*] = [X°] T [C] [X] (2.4) 

* * * 

Where [M ], [K ] and [C ] are the transformed property matrices. 
[M ] and [K ] are diagonal because the generated load-dependent 
Ritz vectors [X°] are [M] -orthonormal and [K] -orthogonal. If 

proportional damping is used. 
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[C] = a [M] + a [KJ (2.5) 

° 1 

where a^ and are constants, then [C*J Is also diagonal. 
Thus, the final system equations of motion are diagonalized and 
reduced to L uncoupled, linear, second order, ordinary 
differential equations, as follows: 

[M*] {2} + [C*l {Z} + [K] {Z} = [X°] T {f (s)> g(t) (2.6) 

These uncoupled equations are solved using any numerical 
technique, obtaining the desired response. Wilson et al. [12] 
also introduced an expression that provides an error estimation 
for the dynamic analysis. As mentioned before, the numerical 
examples introduced by Wilson indicate that by using a fewer 
number of load-dependent Ritz vectors, more accurate results can 
be obtained than by using eigenvectors. 

In a paper by Nour-Omid and Clough [13] in 1984, it was shown 
that by orthogonal! zing any obtained load-dependent Ritz vector 
with only the two preceding vectors (step 3 in Wilson’s 
algorithm), one can theoretically ensure orthogonality with all 
previously calculated vectors. The main problem in this procedure 
is that orthogonality can be lost due to roundoff errors, 
accumulating from one step to another. Thus, orthogonality must 
be checked, while generating each vector, by a certain scheme 
given in reference [13]. Also, it was shown how the loss of 
orthogonality can be corrected if it exceeds a certain tolerance 
limit. The method of generating the load-dependent Ritz vectors 
in reference [13], followed that outlined in Wilson’s algorithm 
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through step number three. Hence, no orthogonal izat ion was 
carried out with respect to the stiffness matrix, hence avoiding 
the solution of the (L x L) eigenvalue problem. The transformed 
stiffness and damping matrices thus obtained are not diagonal. 
However, it was shown in reference [13], that the equation of 
motion (2.2) can be reduced to a tridiagonal form, through a 
prescribed set of transformations. (Note that proportional 
damping, given by equation (2.5), is used). The final transformed 
equation of motion will be an follows: 


[TJ {2} ♦ 


where 


[T 


IX J 


{e> 


[ T l ] + [IJ j {2} + {2} = {e> g(t) (2.7) 

= Transformed mans matrix of order L x L. 

It is a tridiagonal matrix. 

= Transformed stiffness matrix. 

= Identity matrix of order L x L. 

= Transformation of (f(s)} vector. All its 
elements are zeros except the first 
element is non-zero. 


a , a , {2>, {2>, (2> and g(t) ane as before. 

As seen in equation (2.7), the excitation is applied only in the 

first of these equations of motion. While equations 2, 3 up 

to L ane equations of free vibration. 

In another paper by Arnold et al. [14] in 1985, it was stated 
that Wilson’s algorithm had been implemented into the MSC/NASTRAN 
program. The method presented by Wilson et al. [12] was applied 
to structures having 1000 degrees of freedom. The results 
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presented in reference [14] showed that by using load-dependent 
Ritz vectors, instead of eigenvectors, the following was obtained: 

1- A saving of 90 % for eigenvalue extraction. 

2- A saving of 50 % for response calculations. 

In 1986, Wilson and Bayo [15] introduced load-dependent Ritz 
vectors to fixed-interface component mode synthesis, in the 
calculation of forced responses by dynamic substructuring. The 
presented method had the following limitations: 

1- The method is only suitable for forced vibration problems where 
there are external loads, thus it cannot be applied to free 
vibration problems. 

2- The presented formulation is only applicable to fixed- interface 
components. 

3- The displacement shapes used for components are constraint 
modes [4M SLnd load-dependent Ritz vectors for components 
having external loads. The load-dependent Ritz vectors are 
generated from those external loads, assuming the interfaces of 
components totally fixed. Thus, if there are no external loads 
or they are applied at the interfaces of a component, then no 
load-dependent Ritz vectors are generated for this component. 
Hence, only constraint modes will be used, as displacement 
shapes for such a component, and the component representation 
is not expected to be good in that case. Wilson and Bayo 
su gS es t» for such components, to generate load-dependent Ritz 
vectors from a fictitious uniform external load, applied to the 
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component. In the presented problems, the components were 
planar beams and the fictitious uniform load was taken 
perpendicular to the length of the beams. The assumption of 
the uniform load does not have any theoretical basis and cannot 
be justified for complex three-dimensional components, where an 
assumption of a suitable uniform loading is not possible. 




Chapter 3 


Boundary Flexibility Vector Method 
of Component Mode Synthesis 

3 , 1) Introduction: 

The boundary flexibility vector method of CMS is presented in 
this chapter. All of the underlying equations will be derived. 
The method is described for both fixed and free-interface 
components. It should be noted that the developed formulation has 
many similarities to that of conventional CMS methods, thus 
facilitating the combination of conventional normal modes and 
boundary flexibility vectors, as a set of generalized coordinates 
for components. The boundary flexibility Ritz vectors are 
generated by an extension of Wilson’s load-dependent Ritz vector 
algorithm, (described in the previous chapter). The boundary 
flexibility Ritz vectors are not, however, generated from external 
loads, applied to the system. Thus, they can be incorporated into 
both forced and free vibration problems. The generation of the 
boundary flexibility Ritz vectors does not require the solution of 
a costly eigenvalue problem, associated with the equation of 
motion of a component. Hence, a substantial saving in CPU time is 
generally achieved by using the boundary flexibility vectors 
instead of the conventional eigenvectors. 

To illustrate the boundary flexibility vector method of CMS, 
a generic type of structural component will be used, as shown in 
Figure (1). Assume the total number of D.O. F. in the component is 
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n. Let {u } and {u } be the subvectors of internal and interface 

l C 

degrees of freedom in component. The lengths of {u > and {u > are 

i C 

i and c respectively. Hence, the total displacement vector {u} is 
partitioned into 


{ u} 



^ i ' 

Accordingly, the mass [m] and stiffness 
component, can be partitioned into: 


(3. 1) 

[k] matrices of undamped 


[m] = 

[m ] 

CC 

[0 c! ' 

and [k] = 

[k ] 

CC 

Ik J ' 


■ l0 il 

‘",1 . 


[ k 1 

L 1 c 

[ \l - 


(3.2) 


The method will be presented in the following sections, for both 
fixed and free-interface components. 

3 .2) Fixed-Interface Boundary Flexibility Vector Method of CMS: 

As mentioned before, the formulation can be based entirely on 
the boundary flexibility Ritz vectors, or on a hybrid boundary 
flexibility vector/convent ional CMS formulation. The two 
formulations are presented in the following: 

3 .2.1) Boundary-Flexibi 1 i ty Vector Formulation : 

Consider the equation of motion for the undamped 
fixed-interface component 

[mj + [kj {uj = { f (s ) } g(t) (3.3) 

{u^} = acceleration vector of internal D.O.F. 

{ f ( s ) } and g(t) are as before in the previous 
chapter. 


where 
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The reason the equation is written in terms of the internal 
D.O.F. {u } only, is that the interface D.O.F. {u } are assumed 

1 c 

totally fixed, in the case of fixed- interface components. As the 

interface D.O.F. are not actually fixed and they undergo 

displacements within the motion of the whole system, then the 

forcing function on the component can be considered as support 

motions (accelerations) at the interface degrees of freedom. The 

idea of the pseudostatic influence vectors, commonly used in 

earthquake-response analysis (see Clough and Penzien [16]), is 

employed to obtain the loading function. According to reference 

[16], a pseudostatic influence vector defines the static response 

of the internal D.O.F. due to a unit displacement (motion) applied 

at one support (interface) D.O.F. In case each interface D.O.F. 

is subjected to a unit displacement, the total static response of 

the component is obtained by superposition of the pseudostatic 

influence vectors obtained from each independent unit 

displacement. In case the load is defined by a unit acceleration 

at one interface D.O.F., the pseudostatic influence vector is used 

to define the distribution of the load (acceleration) within the 

component. The definition of the pseudostatic influence vectors 

is thus similar to that of the constraint modes [$ ], previously 

c 

defined in chapter (2). The size of matrix [$ ] is (n x c) and it 

c 

contains all the constraint modes in its columns, where j — 

c 

If 2, c. [$ c ] is given by: 
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[* ] 1 


[I 1 1 

= 

cc 

s 

cc 


[*, 1 

L ic J 


[*,] 

u Ic J 


where [I ] = identity matrix of order c 

[♦jl = part of constraint modes defining displacements 
of interior D.O.F. 

[$j] is obtained from solving the lower partitions of the 
equations 

[k ] [k 1 1 

cc cl 

[k] CkJ 

where = interface forces between components 

thus [$ ] = - [k J" 1 [k ] 

ic ii ic 

For an acceleration u^t) given at interface D.O.F. number j, the 
forcing vector will be 

{f (s)> g (t) = - [m l {<p\} ii (t) (3.6) 

J J II 1C CJ 

Thus the total forcing vectors of equation (3.3) is given by: 

{f(s)> g(t) * - £ ( [mj {^} ujt) ) (3.7) 

The summation is carried out for all the accelerations applied at 
all the D.O.F. contained in {u^}. It is clear that for the j th 
forcing vector given by equation (3.6) 

gj(t) = - a c jt) 

Thus the j th spatial distribution of the force is 


[I ] 


[ p 11 

cc 

— 

cc 

[*,] 


. [0.1 . 


(3.8) 
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<f (s)} = [m l {^> (3.9) 

Wilson’s algorithm can then be used to generate the boundary 
flexibility Ritz vectors, from each of the derived spatial 
distributions of forces. The j th starting vector {q*^ is obtained 
from (f^(s)> by 

Ik i! { <S * [ “,1 <*;> (3-101 

Wilson’s algorithm is used without the last step (step 4) of 
orthogonal izat ion with respect to the [k^] matrix. Each 
subsequent vector in the j th set is then orthogonal i zed with 
respect to the mass matrix, with all the previous vectors 
generated in all the sets from 1 through j (step 3 in Wilson’s 
algorithm). Each generated vector is normalized with respect to 
the mass matrix. 

Let l be the number of generated boundary flexibility Ritz 
vectors. Let [Q^] be the matrix containing the generated vectors 
in its columns. The size of [Q £ ] is (i x l ) . Let {p^} be the 
vector containing l generalized coordinates associated with the 
generated boundary flexibility vectors. Then [$ ] and [QJ are 

c c 

used to transform the coordinates, from a physical coordinate 
system to a mixed physical and generalized one. The 
transformation relation is as follows: 


r u] 


f u ] 

r 

Uj 

► = [T3 H 

1 c 1 
l J 

► z 


[I ] 

cc 




[*,] [QJ 

lc l 


(3. 11) 
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where [T] = transformation matrix of order (n) x (c+£) 

It should be noted that the number of boundary flexibility 
vectors, l, must be less than the total number of internal D.O.F. , 
i, in the component, for any reduction of the total number of 
D.O.F. , to be achieved. This is usually the case, an a very small 
number of boundary flexibility vectors, is often sufficient to 
obtain a good representation of a component, as will be 
illustrated in the numerical examples. The transformation matrix 
[T] is used to transform the property matrices, given in equation 
(3.2), to the reduced mixed coordinate system. The final 
transformed component mass matrix [p] and component stiffness 
matrix [Q] are as follows: 





r in ] 

Im], 



(Ml = 

: m T 

[m] [T] = 

cc 

CC 

(3. 

12) 





l “'u ■ 



where [u ] = 

(m ] + 

[*, ] T [m ' 

1 [* J 




cc 

CC 

lc ll 

[ ic 





-P 

O- — 1 
H 

II 

*/ [m i! 

t* ] 

ic 




‘"’a" 

[Q / 

M t<y 

* lli u 

- identity 

matrix 

of 


order l. 


t r ta 1 [q] > i 

[fll = [T] T [k] [T] = cc cl (3.13) 


where [fl ] = [k ] + [k J [• ] * [k ] - [k ] [k J ' 1 (k ] 

CC cc cl ic cc cl II lc 
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IQ.] = cmi = [OJ 

cc c c cc 


IV l \l IV 


For every component, the boundary flexibility vectors are formed, 
and the transformation of coordinates is carried out to obtain 
reduced property matrices. As all displacements of the interfaces 
of the components are defined by physical D.O.F. , then the 
assembly of system property matrices, at the interfaces, is simply 
done by direct summation of the interface portions of property 
matrices of components, as in conventional direct stiffness 
assembly of finite element property matrices. Any type of dynamic 
analysis is carried out using the reduced assembled system 
property matrices. Any obtained results are back transformed to 
the physical coordinate system, by utilizing equation (3.11). 

3 .2.2) Hybrid Boundary Flexibility Vector/ Conventional CMS 
Formulation: 

The method can be applied for the general case where a 
combination of boundary flexibility vectors and natural vibration 
modes (eigenvectors), sire used as displacement shapes of 
components. This case could arise when the number of available 
natural modes of vibrations, obtained from experiments, for 
example, are insufficient to give adequate representation of a 
component. 

Let l be the number of obtained boundary flexibility vectors, 
k be the number of any kept natural vibration modes, [fc ] be the 




26 


matrix containing the kept natural modes in its columns, and {p^} 
be the vector of generalized coordinates associated with the 
natural modes. [Q^l and {p^> are as before. It should be noted 
that the generated boundary flexibility vectors, contained in 
[Q^], are orthogonal to each other and they are orthogonal i zed 
with all the natural modes contained in [♦ ], with respect to the 
mass matrix [m l. [Q.l, [<f ] and [$ ] are used to transform the 

i I l ic k 

coordinates from the physical coordinate system to the reduced 
mixed physical and generalized one. The transformation is as 
fol lows: 


f u] 

, 

r 

U 

C 


ti i 

[0 h 

[01 1 


f 

u 

c 

0 

• = [T 1 i 

p * 

L p k - 

► — 

cc 

c k 


Pp 

UJ 

1 


[• i 

L ic 

[Q,l 

1*1 
k J 


c 

p k J 


where [T ] = transformation matrix of order (n) x (c+i+k) 

The property matrices of equation (3.2) are transformed, giving 
the following reduced mass [/x] and stiffness [Q] matrices. 


[JJ] = [tV [m] [t'] 




[p cl 1 


l » ] u 

tfX d 

o. — . 

Hi 

Kl J 


(3. 15) 


where [/lx ] = [m ] + [* T ] [ml [$ ] 

cc cc ic if ic 

- '"i; * >v T [ ",i 

‘“J - ig T = ig T [»| [« i 
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l *u * "V T [ *,l 'V 

= [1]^ = identity matrix of order l 

'“le ‘ ll *el T ‘ lo le 

l \l = 'V T >",! [ V 

= [I kJ | = identity matrix of order k 


[Q] = [t'] T [k] 


where [£2 ] = [k ] + 

cc cc 


[£2 ] 
CC 

in h 

“U ' 


[a (i 


[a d 

(3. 16) 

la] 

*■ kc 


“U J 


1 - [k ] 

c CC 

- tk cl 

1 [k ,r 

1 [k i ] 
ic 


[s y - > n ]5 - i°(j 


to I - [n 1 T - [01 

kc ck kc 


[nl a - IQ ( ] t (k ] [Q i J 

[a Je ' ‘ a <T “ !*,! I0<1 

£Q kl | = t$ k l T t* k ] = diagonal matrix, having the 


squares of natural frequencies in its diagonal. 
The rest of the procedure is similar to what was explained in the 
previous section. 

3 .3) Free- Interface Boundary Flexibility Vector Method of CMS : 

Two formulations are presented; formulation based only on 
free- interface boundary flexibility vectors and a hybrid 
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free-interface boundary flexibility vector/conventional CMS 
formulation. 

3 .3.1) Boundary Flexibility Vector Formulation : 

Analogous to conventional method of free-interface CMS, two 
cases have to be considered separately; systems in which 
components have rigid-body modes and systems in which components 
have no rigid-body modes. 

3 .3. 1.1) Components Having No Rigid-Body Modes : 

The equation of motion for a n undamped free-interface 
component in physical coordinates is 


1 

1 — 1 
3 

i— i 

[0 ] 1 

aw 


f 5 .l 

► + 

1 

[k ] 1 

aw 

•< 

r y 

u 

a 

N- S . 

f f .l 

| 

[0 ] 
u wa 

[m ] 

WW J 

' 



[k ] 

L wa 

[k ] 

WW J 


u 


L o . 

V. w J 


where {u^} = subset of total displacement vector {u} for which 
internal or external forces are applied. 

{u > ■ complement of {u > in {u>. 

W a 

{f^} = vector of external and internal forces acting on 
component. 

Without a loss of generality, the case of free vibration of a 
component will be considered. The method that will be presented 
in this section can also be applied to the case of components 
having external loads, as will be explained in another section. 
In the case of free vibration, there are only internal forces 
acting at the interfaces between components. Thus {u } is equal 
to {u^} and equation (3.18) will be 
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(3. 18) 



where {f^> = subvector of internal forces acting on the 
interfaces of component. 

The static response of the component to any individual force (or 
moment) in {f } acting at one of the degrees of freedom of the 

C 

interface {u >, is a multiple of its response to a unit force (or 

C 

moment) applied at this degree of freedom. According to the 

previous chapter, the static responses of the free- interface 

component to a set of unit forces (or moments) applied 

individually and successively on every D.O.F. in {u }, were 

c 

defined as the attachment modes [$ ] . The number of attachment 

a 

modes in the component are equal to the number of the degrees of 
freedom in {u^}, where the unit forces are applied. Hence the 
dimension of [$ ] is (n x c) and it contains all the attachment 

a 

modes {^> in its columns, where j = 1, 2, c. The static 

response of the component to the vector {f > is a linear 

c 

combination of the attachment modes in [4 ]. According to the 
definition of attachment modes, they are the columns of the 
flexibility matrix ( [g] = [k]" 1 ) corresponding to the D.O.F. in 
{u } (or {u } in case of free vibration). Hence 

a c 


[* 1 


[4 ] 1 

a cc 

[4 ] 

a ic J 



(3. 19) 
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The free- interface boundary flexibility Ritz vectors are 
derived from the force vector {f }. They are obtained from the 

C 

displacement shapes of the static responses of the component to 
the forces (or moments) in {f>. The j th starting boundary 

C 

flexibility vector, (q^ is obtained from the inertial loading of 
the static response {<f>h , as follows: 

a 

[k] {q*> = [m] (<£ J > (3.20) 

l J a 

Note that the j th starting vector is not taken as {^ J >, which is 

a 

the static response to the applied unit load at interface D.O.F. 
number j. This approach is different than Wilson’s method, where 
the starting vector is taken as the static response to the 
external applied load. This approach is also different than that 
presented in fixed- interface boundary flexibility method of CMS, 
where the starting vector was taken a s the static response to the 
load derived in equation (3.9). The reasons for this approach 
will be clarified in a subsequent section of the thesis. 

After obtaining the j th starting vector, subsequent vectors 
are also obtained in a somewhat different procedure than that of 
the fixed- interface algorithm. The difference is that all the 
normalizations and orthogonal izat ions of vectors are performed 
with respect to the stiffness matrix instead of the mass matrix. 
The reason for this difference will also be clarified in a 
subsequent section of the thesis. It should be noted that any 
obtained vector In the j th set is orthogonal ized with respect to 
the stiffness matrix, with all previous vectors obtained in sets 1 
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through j. 

In a paper presented by Hintz [17], it was pointed out that 
in order to obtain a good representation of components in 
free-interface CMS, an accurate determination of displacements at 
D.O.F. in {u } (or {u > in free vibration), where forces are 

a c 

applied, is required. In other words, the static displacement 
response of the component to interface forces must be the same for 
both the CMS model an d an acceptable finite-element model of the 
component. In this case, the static representation of the 
interface D.O.F. is said to be complete, in the new system of 
generalized coordinates. In order to obtain the same displacement 
response, for both the CMS and finite-element models, the 
representation of the columns of the flexibility matrix 
corresponding to the D.O.F. in {u^} must be complete for the CMS 
model. If the complete set of free-interface normal modes is used 
as displacement shapes for the component, then the representation 
of the flexibility of interface D.O.F., is by definition, 
complete. Since a truncated set of the free-interface normal 
modes is generally used in conventional CMS, the flexibility of 
the component is generally less than that of the corresponding 
finite-element model. That is, the contribution of the truncated 
set of free-interface normal modes to the flexibility matrix is 
incomplete, and the free-interface normal modes need to be 
supplemented with static displacement shapes. The static 
displacement shapes used in conventional CMS are the residual 
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flexibility modes [$ d ], defined before in chapter two (see 
references [2], [5], [6], and [17]). The residual flexibility 

modes [# ] were introduced by MacNeal [5] to replace attachment 
modes [$ ], as static displacement shapes for components. The 
reason for introducing the residual flexibility modes is that by 
using attachment modes [# ] to supplement the truncated set of 
free-interface normal modes, the obtained flexibility 
representation of the interface D.O.F. in the CMS model is greater 


than that of 

the 

finite-element 

model. 

The more flexible 

represent at i on 

is 

attributed to 

the 

complete 

f lexibi 1 ity 

representation 

of 

the interface D. 

O.F. , 

provided 

by using the 


attachment modes as the only displacement shapes for components. 
This fact explains why the static responses to the interface 
forces (which are the attachment modes), were not used as the 
starting boundary flexibility Ritz vectors. In the event they 
were utilized as starting vectors, the obtained interface 
flexibility of a component would actually be greater than that of 
the corresponding finite-element model. 

In the free-interface boundary flexibility method of CMS, the 
contribution of the obtained boundary flexibility vectors to the 
interface flexibility is also incomplete. To obtain a complete 
flexibility, the contribution of the free-interface boundary 
flexibility vectors to the flexibility matrix must be determined. 
Assume the number of the obtained boundary flexibility vectors is 
l. Let [Q^] be the matrix containing the l boundary flexibility 
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vectors in its columns and let {p^> be the vector of generalized 
coordinates associated with the free— interface boundary 
flexibility vectors. The sizes of [Q^] and {p^} are (n x £) and 
(£) respectively. Consider the static equilibrium equation of the 
component 


where [k] 
{ u> 


tk] {u} 


f 

i 

0 


v l " 

total stiffness matrix of component, 
total displacement vector. 


(3.21) 


>■ = force 


vector containing only forces 


at 


interface degrees of freedom. 

The boundary flexibility vectors [Q^] are used to transform the 
physical coordinates {u} to the reduced set of generalized 
coordinates {p^}, according to the following 


(u> = [Q £ ] {p £ } 

Substituting equation (3.22) in 
pre-mult iplying by then 

[Q £ ] T tk] [Q t ) {p^ = [<y T 

From which 


equation 


f 1 


°i J 


r r 


-1 

ff 1 

I [Q £ ] tk] 

[Q,J 

] “V • 

O -H 

O 


(3.22) 
(3.21) and 


(3.23) 


(3.24) 


Substituting equation (3.24) in equation (3.22) 
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{u> = [Qj] 


[<y T Ik] [<y 


-1 


“ [q/ 


The static response, {u> , to the applied force vector 


f c 1 

1 


- c 

Y (3.25) 

0 J 

1 J 


f f ] 

force 

vector * 

1 c 

► 



1 0 


is equal to 



(3.26) 


i ' 

where [g^] = flexibility matrix of component represented by 
boundary flexibility vectors. 

Hence by comparing equation (3.25) to equation (3.26), the 
contribution of the free- interface boundary flexibility vectors to 
the flexibility matrix is 


[g k ] = [( V [ [Q £ ]T tk] [( V ] tQ £ ]T (3.27) 

For free-interface boundary flexibility vectors which are 
orthogonal i zed and normalized with respect to the stiffness 
matrix, [g^] becomes 

[gj = [Qj] [Q £ ] T (3.28) 

This relationship explains why the generated boundary flexibility 
vectors were orthogonal ized and normalized with respect to the 
stiffness matrix [k] instead of the mass matrix [m]. Notice that 
the indicated inversion of the matrix in equation (3.27) is 
avoided by the orthogonal izat ion and normalization with respect to 
the [k] matrix. The unrepresented flexibility, or residual 
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flexibility [g ], analogous to that of conventional free-interface 

d 

CMS, is given by 

[gj = fg] ~ [g.l (3.29) 

d k 

where [g] = full flexibility matrix = [k] _1 

[g^] is given by equation (3.28) 

The residual flexibility matrix [g ] is used to obtain the 

d 

residual attachment modes (residual flexibility modes) [¥ ], by 

d 

applying unit forces or moments at D.O.F. in {u > (or (u > in free 

a c 

vibration). Hence, 


[* 1 

d 


' [i ] ' 

[g.l 

[o j 

u ic J 


‘^cc 

[ 6 d l lc 


(3.30) 


where (I ] = identity matrix of order c. 

CC 

= matrix of unit forces or moments applied at 
D. 0. F. in {u > . 

c 

Notice that the residual attachment modes ] are the columns of 

d 

the residual flexibility matrix [g ], corresponding to the degrees 

d 

of freedom in {u }. 

c 

The boundary flexibility vectors [Q^] and the residual 
attachment modes [¥ ], will be used to transform the equation of 

d 

motion from the physical system of coordinates to the reduced 

generalized one. (QJ and [¥ ] will provide the complete 

-c d 

flexibility for the D.O.F. in {u } (or {u } in free vibration), 

a c 

where forces a re applied on the component. Let {p > and {p.} be 

d L 

the generalized coordinates associated with [tf 1 ] and [QJ 

d c 
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respectively. The sizes of {p > and {p.} are c and i 

d t 

respectively. The transformation of coordinates is given by 


K! 

1 

► ■ [T 1 - 


^ = 

[®J 

d cc 

[Q 1, ' 

C t 


[ P <1 

j 1 

(3.31) 
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i J 

1 

l P i J 


[*,1, 

L die 

Jo 

i 


Ip* J 



where [T^J = transformation matrix of order (n) x (c+£). 

[T^ can be used to transform the property matrices of component. 
However coupling of components, whose interface displacements are 
expressed by a set of generalized coordinates and mode shapes, is 
not straightforward as in fixed- interface boundary flexibility 
method of CMS. To overcome this problem, the displacements of the 
interface D.O.F. are back transformed to physical D.O.F. (see 
Abdallah [2]), in order to allow for a direct stiffness assembly 
process. Assume the total number of generalized coordinates is J, 
and if c is the number of interface D.O.F., and m - j - c, then 
let 

= vector of all generalized coordinates = 

(p > = any subset of {p > having a size of c. 
c J 

{p } = complement of {p } in {p }. 

m c J 

Accordingly, the transformation matrix [ T^ ] in equation (3.31), is 
partitioned and rewritten as follows 


/• * 
u 

c 

- » [T 1 - 

r p o i 

_ 

[T 1 

1 CC 

[T ] 1 

1 cm 

•< 

P c 

(3.32) 
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L 1 J 
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LpJ 


[TJ. 

u 1 lc 

IT 3 

1 in J 
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(3.32) 
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From the first matrix equation in (3.32) 

{p > = IT I" 1 {u > - [T ] _1 (T ] {p > 

c X cc c 1 cc 1 cm m 

And from the second matrix equation in (3.32) 


(3.33) 


<u > = [T ] {p > + [T ] {p > (3.34) 

l 1 lc C 1 ll ■ 

Substitute equation (3.33) in equation (3.34), then 

<u } = [T 1, [T l~ l <u > + ( [T ] - [T ] [T I" 1 [T 1 ) <p > 

l 1 lc 1 cc c 1 !■ 1 lc 1 cc 1 cm m 


(3.35) 

By using equation (3.35), the transformation of coordinates in its 
final form is as follows 


K1 

■ = [T ] • 
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li 1 

cc 
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( u =i 

K J 
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L p J 

w m J 
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[#*] 
L ic 
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i m J 
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(3.36) 


where [T^I = transformation matrix of size (n) x (j) 
[**] = [T ), (T, I" 1 

lc 1 ic 1 cc 


[QJ = [T , ] t * ITJ, f T J _1 tTl 

lm 1 In 1 lc 1 cc 1 cm 


Notice that the form of equation (3.36) is similar to that of 
equation (3.11), for f ixed-interface components. The property 
matrices are transformed to the reduced coordinate system by using 
the transformation matrix (T ] , of equation (3.36). As the 
interface D.O.F. are defined with respect to the physical 
coordinate system, then the assembly of system property matrices 
is straightforward and similar to that performed for 
f ixed-interface components, in the previous section. 
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3 .3. 1.2) Components Having Rigid-Body Modes : 

The main complication of components having rigid-body modes 
is that the stiffness matrix is singular; thus the flexibility 
matrix cannot be directly obtained. This problem was Investigated 
by Rubin [6] in conventional CMS. He introduced an alternative 
method for obtaining the residual attachment modes. To overcome 
the singularity of the component stiffness matrix, he used inertia 
relief loading (see references [2], [4], [6] and [7] for a 
detailed explanation of the method). The results obtained by 
Rubin [6] are stated briefly here. 

Assume the number of rigid-body modes is r. Let [4 ] be the 

r 

matrix containing the r rigid-body modes in its columns. The 
rigid-body modes are normalized with respect to the mass matrix. 
Let {u} be the vector of all physical D.O.F. in component and let 
{ u } be a subset of {u} where forces are to be applied to obtain 
attachment modes. In Rubin's method, the component is constrained 
and prevented from rigid-body motion by applying r constraints to 
a set {u > of D.O.F. , such that {u > is exclusive of {u }. The 

r r a 

sizes of {u}, {u } and {u } are n, a and r respectively. Let {u > 

a r w 

be the complement of {u } and {u } in {u>, having a size of w. 

a r 

Define a special flexibility matrix [g 3 relative to {u > for the 

c r 

constrained component. 
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Define a square projection matrix [A] as follows 

[A] = [I ] - [m] [« ] ] T 

nn r r 

where [I ] = identity matrix of order n. 

nn 

Rubin showed that the elastic flexibility matrix [g ] is given by 

[g 1 = [A] T [g ] [A] (3.38) 

e c 

The flexibility matrix [g^] is used to obtain inertia relief 
attachment modes [¥ ], by applying unit forces at {u } (or {u > in 

a a c 

case of free vibration). Inertia relief attachment modes are the 
columns of the flexibility matrix [g ] corresponding to the 
degrees of freedom in {u }. 

a 

Similar to components having no rigid-body modes, where 
attachment modes [$ ] were used to generate the boundary 
flexibility vectors, the Inertia relief attachment modes [¥ ] are 

a 

used here to generate the starting vectors. The j th boundary 
flexibility starting vector {q*>^ is obtained from 

(q*> , = tg 1 Cm] (^ J > (3.39) 

1 J a a 

where = j th inertia relief attachment mode. 

[g^l = flexibility matrix given by equation (3.38). 

The subsequent boundary flexibility vectors are normalized and 
orthogonal! zed with respect to the stiffness matrix, with all 
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preceding vectors, as in the previous section. The contribution 

of the boundary flexibility vectors to the flexibility matrix is 

as before and given by equation (3.28). Hence the residual 

flexibility matrix [g 1 is obtained from 

a 

[g d ] = [g e ] ” [( V [Q * ]T (3.40) 

The residual flexibility matrix [g ] is used to obtain residual 

d 

inertia relief attachment modes [¥ ], by applying unit forces at 

S 

{u > (or {u } in case of free vibration). The residual inertia 

a C 

relief attachment modes are the columns of the residual 
flexibility matrix [g ], corresponding to the D. 0. F. in {u >. 

d a 

($ r ), [Q^] and [4^] are used to transform the physical coordinate 
system to a generalized one as follows: 


/■ 

u 

a 


t* ] 

as 

[Q], 

[* 1 1 
a r 


p 

P s 

- U 

- = 

r* ) 


[• ) 


Pi 

w 


W3 

wr 



U 

r J 


] 

L rs 

tQ h 

i* ] 

r r J 


IK J 


where {p > = generalized coordinates of [¥ ], of size a. 

s s 

= generalized coordinates of [Q^], of size Z. 

(p } = generalized coordinates of [$ ], of size r. 

r r 

The size of the transformation matrix of equation (3.41) is 
(n) x (a + l + r). The transformation of property matrices and 
the assembly of the system property matrices will be similar to 
what was presented in the previous section. 
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3 .3.2) Hybrid Boundary Flexibility Vector/ Conventional CMS 
Formulation : 

In this case a truncated set of free- interface normal modes 
[4 fe ] is used to supplement the free-interface boundary flexibility 
vectors, as displacement shapes for components. Let {p } be the 

k 

vector of generalized coordinates associated with [$ ]. The sizes 

k 

of {p^) and are (k) and (n x k) respectively. The following 

steps are general regardless the component has rigid-body motion 
or not. 

1- Normalize all the free-interface natural modes [$ ] with 

k 

respect to the component stiffness matrix [k]. ( optional step). 

2- Generate the free-interface boundary flexibility vectors [Q^] 
from the starting vectors of equation (3.20), in case of 
components having no rigid-body modes, or equation (3.39), for 
components having rigid-body modes. It should be noted that 
every generated vector is orthogonal ized with respect to the 
stiffness matrix [k], with all natural modes contained in [$ ] 

k 

and with all previously generated boundary flexibility vectors. 
Every generated vector is normalized with respect to the 
stiffness matrix. Let t be the number of the generated 
boundary flexibility vectors and let {p^} be the vector of 
generalized coordinates associated with them. 

3- Form a matrix [V] containing the normal modes of vibration [$ ] 

k 

and the generated boundary flexibility vectors [Q ], in its 
columns. The size of [V] is (n) x (k + l) . Hence, 
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[V] = [ [* k I [Q^ ] (3.42) 

Accordingly the vector of total generalized coordinates, {p }, is 

V 

formed as follows 


{p v } 



(3.43) 


The contribution of the [V] matrix to the flexibility matrix [gj 
is 


(gj = [V] [V] T (3.44) 

Or in case the natural modes of vibration [$ ] are not normalized 

k 

with respect to the [k] matrix (step 1 above), [g^] is given by 
(see equation (3.27)) 


r t v 1 T 

[gj = [V] [V] T [kj [V] [V] T = [V] 



[0 h 

■ I0 d 

[l] u 


1-1 


[V] T 


(3.45) 

where [A ^ = a diagonal square matrix of order k. It 

contains the eigenvalues, corresponding to 
natural modes [• ] f in its diagonal. 

[I * identity matrix of order L 
Although equation (3.45) involves an inversion of a matrix, it can 
be used because the matrix is diagonal and its inverse can be 
obtained readily. 

The subsequent steps, whether for components with or without 
rigid-body modes, are similar to that explained in the previous 
section, except that every [Q^] is replaced by [V] and every {p^} 
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is replaced by {p^}, in the equations. 

3 .4) Computational Effort for Boundary Flexibility Vectors vs 
Normal Modes of Vibration : 

In this section, the number of operations required to 
generate U) boundary flexibility vectors versus (k) normal modes 
of vibration, is estimated. In the numerical problems of chapters 
four and five, the method used for obtaining the eigenvalues is 
the determinant search method (see Bathe [18]) with bisection 
iteration. While the method used to extract the mode shapes 
(eigenvectors) is the inverse iteration method (see Bathe [18]). 
The number of required operations is estimated by using some basic 
criteria from reference [19]. The detailed calculations are 
listed in appendix A, while in this section only the final results 
are presented here. 

It was found that the number of operations required to obtain 
(k) eigenpairs, for a component, is as follows: 

Total number of multiplications and divisions = 

(3.46) 

Total number of additions and subtractions = 

-jp (8n 3 -9n 2 +n) + £ £ -p - (n 3 +3n 2 -n) + x t (2n 2 +n- 3) J (3.47) 

where n = total number of physical D.O.F. in component. 

= number of inverse iterations required to obtain 
eigenvector i. 


k 

3 


(4n -n) + 


k r y 

i v 

i = 1 L 


(2n 3 +9n 2 -5n) + x t (2n 2 +4n +1) 
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Yj = number of iterations, in determinant search 
method, required to obtain eigenvalue i. 

While for the generation of l boundary flexibility vectors per 
component, the number of operations is as follows: 


Total number of multiplications and divisions = 


x 3 A , l 2 + 51 , 2 

) n + ( — = ) n + 


, 31 +31-1 . 

( = ) n (3.48) 


Total number of additions and subtractions = 



n 3 * ( ^ ?* - 3 ) n 2 + ( 3 ' 2 - 13 ' 
d o 

+ 1 1 
) n 

. l 2 +t 
1 2 

) (3. 

49) 

It can be 

seen from 

equations (3.46) 

through 

(3.49), 

that 

the 

number of 

operations 

is proportional to n 3 . 

Thus in case 

of 

components 

having a 

large number of 

D.O.F. 

, the 

number 

of 

operations 

will be 

governed by the 

n 3 terms in 

the above 

equations. 

By using 

only the highest 

order cubic terms in 

the 


above expressions, approximate values for the number of operations 
can be obtained as follows: 

- For the generation of k eigenpairs: 

Approximate number of multiplications and divisions = 

( 4 + y ) ( JL. ) n 3 (3.50) 

where y = average number of iterations in determinant search 


method. 
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The approximate number of additions and subtractions is found to 
be the same as in equation (3.50). 

~ For the generation of l boundary flexibility vectors: 
Approximate number of multiplications and divisions = 



Also the approximate number of additions and subtractions is the 
same as in equation (3.51). 

It is worth noting that the number of generated boundary 
flexibility vectors , i, does not appear in equation (3.51). 
Hence, once the stiffness matrix is inverted to obtain the 
starting vector, the subsequent vectors Eire generated almost 
without any computational cost. 

The ratio of the number of operations to extract k eigenpairs / 
number of operations to generate l boundary flexibility (B. F. ) 
vectors can be approximated by: 

Approximate # of operations for extraction of k eigenpairs _ 

Approximate # of operations for generation of l B. F. vectors ~ 

( 1 + ) k (3.52) 

The approximated ratio of equation (3.52) is dependent only on k 

and y and it can be seen that it is always greater than unity. 
Hence the number of operations required to extract k eigenpairs is 
always greater than the number required to generate l boundary 
flexibility vectors. Furthermore, for more complex components, 
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where one would require a substantial number of normal modes (k) 
for a realistic characterization, the ratio given in equation 
(3.52) will be relatively high. 

3 .5) Summary : 

A new method of CMS namely, the boundary flexibility vector 
method of CMS, was presented in this chapter. The displacement 
modes used in the method, represent a new type of vector, named 
boundary flexibility Rltz vectors. Those vectors are generated 
from internal forces derived from the flexibility (or stiffness) 
properties of the interface degrees of freedom of components. It 
should be pointed out that the presented method avoids the 
solution of a costly component eigenvalue problem, which is 
required in conventional methods of CMS, to obtain the natural 
vibration modes of the components. 

A comparison of the number of operations required to obtain 
natural modes of vibration versus the boundary flexibility Ritz 
vectors was made. The comparison indicated that a substantial 
reduction of the number of required operations was obtained by 
using the boundary flexibility vectors instead of natural modes in 
CMS. 

It should be noted that, in all of the methods of fixed or 
free- interface boundary flexibility method of CMS outlined above, 
for the case where there are external loads acting at internal 
D.O. F. {u t > of components, those forces are used to obtain 
load-dependent Ritz vectors (as in Wilson’s algorithm presented in 
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chapter (2)). For this case, the formulation of the boundary 
flexibility method of CMS will be presented in chapter five. 



Chapter (4) 


Results of Numerical Problems 
in Free Vibration 

4 . 1) Introduction : 

This chapter presents three sample problems of free vibration 
solved numerically by the boundary flexibility Ritz vector method 
of CMS (B.F. Ritz CMS). The problems are also solved by 
conventional CMS and finite-element methods. The results obtained 
by the boundary flexibility method of CMS are compared to the 
other methods. Also a comparison of the saving in CPU time, 
obtained by using the boundary flexibility vectors instead of 
eigenvectors, is presented in this chapter. The first two 
problems are composed of small components, having a few number of 
physical D.O.F. While the third problem is subs t suit i ally larger, 
and demonstrates very well the saving in CPU time obtained by 
using the boundary flexibility method of CMS. 

4 .2) Sample Problem One : 

Sample problem one consists of a beam fixed at both ends and 
supported at mid-span by a spring, discretized as shown in Figure 
(4.1) (all units are lbs and inches). The properties of the beam 
are given in the figure. Each node has two D.O.F. ; vertical 

displacement and flexural rotation. 

The eigen-problem for this problem was solved for eleven 
different values of spring stiffness k, ranging from 0.0 to 
50000.0 lb/in. For each value of k, five mode shapes and natural 
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frequencies were obtained by solving the problem by nine different 
methods. The first solution for the problem was obtained from a 
finite-element model of the system. The problem was then solved 
using CMS by dividing the system into two identical components. 
One component was the beam from node 1 to node 9 and the other 
from node 9 to node 17. The constraint and attachment modes, of 
component number one, used to generate the boundary flexibility 
Ritz vectors are shown in Figures (4.2) and (4.3) respectively. 
The various fixed- interface and free- interface CMS methods, 
utilized in solving the problem, are listed in the following 
sections. 

4 .2.1) Solution by Fixed- interface Methods of CMS : 

Four methods were used to solve the problem: 

Method 1- Conventional CMS, using two fixed— interface normal modes 
and two constraint modes, for each component. 

Method 2- Boundary flexibility CMS, using two fixed-interface 
boundary flexibility Ritz vectors and two constraint 

modes, for each component. 

Method 3- Boundary flexibility CMS, using three fixed-interface 

boundary flexibility Ritz vectors and two constraint 

modes, for each component. 

Method 4- Mixed conventional and boundary flexibility CMS, using 
one fixed- interface normal mode, one fixed-interface 
boundary flexibility Ritz vector and two constraint 

modes, for each component. 
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For all methods, the values of the obtained results for 
natural frequencies were compared to that obtained by the 
finite-element method. The percentage of error was computed for 
all the obtained frequencies, over the given range of spring 
stiffness. The computed percentage errors, for all methods of 
solution, for the first two natural frequencies were less than 
0.06 % and these results are not listed. For the third, forth and 
fifth natural frequencies, the comparative results are shown in 
Figures (4.4), (4.5) and (4.6), respectively. As indicated from 

the figures, by increasing the spring stiffness, the % errors are 
generally increasing, except for mode shapes having a node at the 
spring location. It can be seen from the figures, that the 

results obtained by conventional CMS are better than those 
obtained by boundary flexibility CMS with two boundary flexibility 
Ritz vectors per component. It should be noted, however, that the 
third natural frequency computed by using boundary flexibility 
CMS, with two boundary flexibility Ritz vectors per component, are 
quite reasonable for engineering purposes. Moreover, a 
substantial saving in CPU time was attained by using two boundary 
flexibility Ritz vectors rather than two normal modes, for each 
component, as will be discussed in the next section. A 
substantial reduction in CPU time can even be obtained when using 
three boundary flexibility Ritz vectors per component, compared to 
two normal modes (as will be seen in the next section). The 

results obtained by the mixed conventional/ boundary flexibility 
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CMS were generally in between those obtained by conventional and 
boundary flexibility methods of CMS, with two boundary flexibility 
vectors per component. For certain values of spring stiffness, 
however, they were inferior to the boundary flexibility method 
(for the third and forth natural frequencies) and superior to 

conventional CMS (for the fifth natural frequency). The results 
obtained by boundary flexibility CMS, with three boundary 
flexibility Ritz vectors per component, were in better agreement 
with the finite-element method solution than those obtained by 

conventional CMS. 

4 .2.2) Solution by Free- Interface Methods of CMS : 

Four methods were used to solve the problem: 

Method 1- conventional free-interface CMS, using two 

free-interface normal modes and two residual attachment 
modes per component. 

Method 2- Free-interface boundary flexibility CMS, using two 

free-interface boundary flexibility Ritz vectors and two 
residual attachment modes per component. 

Method 3- Free-interface boundary flexibility CMS, using three 
free-interface boundary flexibility Ritz vectors and two 
residual attachment modes per component. 

Method 4- Mixed conventional and boundary flexibility CMS, using 
one free-interface normal mode, one free-interface 
boundary flexibility Ritz vector and two residual 
attachment modes per component. 
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As before, the percentage errors of the first two natural 
frequencies, obtained by all four methods, are very small (less 
than 0.075 % ) and are not listed. The comparative results for 
the third, forth and fifth natural frequencies are shown in 
Figures (4.7), (4.8) , and (4.9) respectively. Comparing the 

results to that obtained by fixed- interface components, it can be 
seen that for conventional and mixed methods of CMS, the results 
obtained by using fixed-interface components were generally better 
than that obtained by using free-interface components. While for 
the boundary flexibility method of CMS (whether two or three 
boundary flexibility vectors are used per component), the results 
were nearly the same for both fixed and free— interface components. 
It can be seen from Figures (4.7), (4.8) and (4.9) that the 

results obtained by conventional free-interface CMS were better 
than those obtained by the boundary flexibility method, with two 
boundary flexibility vectors per component. However, they were 
comparable and the difference was not as large as in 
fixed- interface methods. The results obtained by the mixed 
conventional/ boundary flexibility model were always in between 
those obtained by conventional and boundary flexibility (with two 
vectors per component) methods of CMS. The results obtained by 
the boundary flexibility CMS, with three boundary flexibility Ritz 
vectors per component, were in the best agreement with the 
finite-element method solution, and they were nearly equal. It 
can be seen from Figure (4.9) that the conventional CMS, boundary 



53 


flexibility CMS (with two vectors per component), and mixed 
methods gave unreasonable values for the fifth natural frequency. 

CPU time was compared for the extraction of normal modes of 
vibration versus the generation of the boundary flexibility 
vectors, for the free- interface solutions. For each component, it 
was found that a saving of 47 Y, of CPU time was achieved in 
generating two free- interface boundary flexibility vectors rather 
than two free-interface normal modes. Moreover a saving of 31 % 
of CPU time was achieved, for the case of generating three 
boundary flexibility vectors versus two normal modes. It should 
be noted that the most costly operation in generating the boundary 
flexibility vectors is the inversion of the full stiffness matrix 
[k], in case of free-interface components (see equation (3.20)), 
or the stiffness matrix [k^ of internal D.O.F. , in case of 
fixed-interface components (see equation (3.10)). This inversion 
of the stiffness matrix, is also required in conventional 
free-interface CMS (see equations (3.29) and (3.30)) for the 
generation of residual attachment modes, or in conventional 
fixed- interface CMS (see equation (3.5)) for the generation of 
constraint modes. As this fact was not taken into consideration 
in computing the saving in CPU time, the above results regarding 
the saving in CPU time can be considered conservative. 

4 .3) Sample Problem Two : 

This problem is used to illustrate the solution of free 
vibration problems by free-interface boundary flexibility vector 
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method of CMS, for components having rigid-body modes. The 
problem is similar to problem one with the support at node one 
being hinged instead of being fixed (see Figure (4.10)). Hence 
component number one has one rigid-body mode, which is the 
rotation of the whole component about the hinged node. For a 
fixed spring stiffness, the problem is solved by four methods. 
The first method of solution is by the finite-element method, 
while the other three solutions are obtained by CMS. The three 
methods of CMS used, are as follows: 

Method 1- Conventional free- interface CMS, using one rigid-body 
mode, one free-interface normal mode and two residual 
inertia relief attachment modes for the first component. 
Two free-interface normal modes and two residual 
attachment modes were used for the second component. 
Method 2- Boundary flexibility CMS, using one rigid-body mode, one 
free-interface boundary flexibility Ritz vector and two 
residual inertia relief attachment modes, for the first 
component. Two free-interface boundary flexibility Ritz 
vectors and two residual attachment modes were used for 
the second component. 

Method 3- Boundary flexibility CMS as in method 2 above, but the 
number of free-interface boundary flexibility vectors, 
in the two components, is increased by one. 

Five natural frequencies were obtained for each method. The 
% discrepancy of all the obtained natural frequencies, compared to 
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the finite-element solution, were computed for all the three 
methods (see table (4.1)). For conventional CMS (method 1), the 
results were slightly better than those obtained by boundary 
flexibility CMS (method 2). Both of the two methods gave 
unreasonable values for the fifth natural frequency. The results 
obtained for all five natural frequencies by boundary flexibility 
CMS (method 3) were in the best agreement with the finite-element 
method solution. 

4 .4) Sample Problem Three : 

This sample problem is used to verify the accuracy of the 
results obtained by the boundary flexibility method of CMS, in 
case of components having a large number of D.O.F. The problem is 
also used to compute the saving in CPU time obtained by using 
boundary flexibility CMS instead of conventional methods. 

The problem consists of a three dimensional three storey 
frame, as shown in Figure (4.11), fixed to the ground at four 
columns. The properties of the beams and columns are given in 
Figure (4.11). Also given is the finite-element model for the 
frame. The total number of physical D.O.F. for this frame is 648, 
corresponding to 6 D.O.F. per node. The problem was first solved 
by the finite-element model obtaining the first twelve natural 
frequencies. It was then solved by the boundary flexibility and 
conventional methods of CMS. The structure was divided into four 
components by two vertical planes; the first plane passing through 
nodes 29, 30, 31, 82, 83 and 84 and the second plane passing 
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through nodes 26, 27, 28, 85, 86 and 87. Component number one is 
shown in Figure (4.12). Fixed- interface components were used for 
the two CMS methods. For each component, 36 constraint modes were 
used, corresponding to six D.O.F. per node at the interfaces. For 
the boundary flexibility method of CMS, three boundary flexibility 
vectors were generated, for each component. The same number of 
eigenvectors were used in conventional CMS, for each component. 
Hence the total number of D.O.F. for the whole frame in CMS 
methods was reduced to 84, compared to 648 D.O.F. in the 
finite-element model. The results obtained for the first twelve 
natural frequencies, by the boundary flexibility and conventional 
methods of CMS were compared to that obtained by the 
finite-element method, and the results were presented in Table 
(4.2). From Table (4.2), it is clear that the results obtained by 
both methods of CMS, for the twelve natural frequencies, are in a 
very good agreement to that obtained by the finite-element method. 

The percentage saving in CPU time obtained by using the 
boundary flexibility vectors instead of eigenvectors, was computed 
for each component. Figure (4.13) shows the percentage saving in 
CPU time, in the case where equal numbers of boundary flexibility 
vectors and eigenvectors were generated, up to eight boundary 
flexibility vectors or eigenvectors per component. Two curves are 
shown in Figure (4.13). The lower curve of Figure (4.13) is the 
saving computed for the case where the CPU time required to invert 
the stiffness matrix in the boundary flexibility method of CMS, is 
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taken into consideration. The upper curve of Figure (4. 13) 
represents the saving when the CPU time required to invert the 
stiffness matrix is not taken into consideration. The reason for 
showing the upper curve is that the stiffness matrix inversion is 
also required in conventional method of CMS to generate the 
constraint modes. An average percentage savings in CPU time of 
75/« in the lower curve and 80% in the upper curve are observed. 
Figure (4. 14) shows the percentage saving in CPU time obtained by 
using up to five boundary flexibility vectors versus one 
eigenvector. For the case where the CPU time required to invert 
the stiffness matrix in the boundary flexibility method is 
included, it can be seen from Figure (4.14) that 12% saving in CPU 
time can be achieved, if four boundary flexibility vectors are 
generated instead of one eigenvector. While for the case where 
the time required to invert the stiffness matrix is not included, 
a saving of 20% in CPU time car be achieved, even if five boundary 
flexibility vectors are generated instead of one eigenvector. 
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Method 

"l 

U 2 

“3 

"4 

3 

1 

0.05853 

0.02541 

0.63544 

2.33334 

25. 0631 

2 

4. 24369 

4.23759 

1.73601 

2.80153 

31.4793 

3 

0.05878 

0.02421 

0.03567 

0.14453 

2. 54257 


Table (4. 1) — A Error of Natural Frequencies For Sample Problem Two 
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Natural 

Frequency 

(rad/sec) 

Method 

FEM 

Conven. 

CMS 

B.F. CMS 

"i 

24.5962 

24.6043 

24.6025 

U 2 

29.6600 

29.6670 

29.6691 

W 3 

30. 2068 

30.2142 

30.2163 

U 4 

30.8600 

30.8678 

30.8698 

U 5 

32. 8940 

32.9237 

32.9275 

W 6 

37. 3674 

37.3866 

37.3686 

*7 

38.6329 

38.6544 

38.6339 

W 8 

39.3779 

39.4619 

39.4755 

W 9 

44.9766 

45.0620 

45.0627 

MEM 

52.0542 

52. 1692 

52. 1414 

0) 

n 

64.4975 

64. 9022 

65. 1329 

0) 

12 

65.8761 


■HH 


Table (4.2) 


Comparison of Natural Frequencies 
for Sample Problem Three 


































Chapter (5) 


Forced Vibration Response by The Boundary 
Flexibility Vector Method of CMS 

5 .1) Introduction : 

This chapter presents the formulation of the boundary 
flexibility method of CMS, for the forced vibration problem. It 
also presents a derivation indicating that the load-dependent Ritz 
vectors, used by Wilson in his formulation for fixed-interface 
CMS, yield a flexibility representation higher them the 
corresponding finite-element model, for the internal D.O.F. , where 
external loads are applied to generate the vectors. It is further 
shown in this chapter, how the load-dependent Ritz vectors can be 
generated in a way such that the flexibility representation of 
fixed- interface components is equal to the corresponding 
finite-element model. Finally, a numerical example is solved by 
several methods and for different loading conditions. The results 
of the boundary flexibility method is compared to that obtained by 
the other methods. 

5 .2) The Application of The Boundary Flexibility Method of CMS in 
Forced Vibration Problems : 

The boundary flexibility method of CMS, presented in chapter 
three, can be applied to forced vibration problems. For 
components having no external loads applied, the formulation will 
be similar to that presented in chapter three, for fixed and 
free-interface components. For components where external loads 
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are applied, the following two cases are considered: 

1) If the external load is applied at an interface D.O.F. , the 
formulation will be exactly the same as to what was presented in 
chapter three, for fixed and free-interface components. 

2) If the external load is applied at an internal D. O.F. , then the 
external load vector is transformed to the mixed physical and 
generalized coordinate system as follows: 

Consider the equation of motion for a fixed or free - interface 
component : 


[m] {u> + [c] {u> + [k] (u> 



(5.1) 


^ l ' 

where {P^} = subvector of interface internal forces between 
components. 

= subvector of external loads applied at one or 
more internal D.O.F. 

[m], [k], [c], {u>, {u} and {u> are as before. 

The equation of motion is transformed to the mixed physical and 
generalized coordinate system as follows: 


[T] T [m][THu> + [T] T [c][THu> + [T] T [k][T]{u> = [T] T < 


f P 1 

c 


1 ' 


(5.2) 

where [T] = transformation matrix given by either of 

equations (3.11), (3.14) or (3.36). 

It should be noted that the internal force subvector, {P >, will 

C 

be unaffected by any of the transformations given by equations 
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(3.11), (3.14) or (3.36). While the external load subvector is 

transformed from the physical coordinate system to a mixed 
physical and generalized coordinate system. 

For free-interface components, the method will be similar to 
that outlined in chapter three, with every {u > vector replaced 
by {u } and every (u } vector replaced by {u >, in all the 

a i w 

equations , where {u }• is the vector of D. 0. F. where external or 

a 

internal loads are applied and {u > is the complement of (u > in 

w a 

{u>. Thus the load-dependent vectors will be contained in the 
residual attachment modes and the boundary flexibility vectors. 

For fixed-interface components, the formulation has to be 
modified in order to generate residual attachment modes from the 
external applied loads. Assume that the external loads are 

applied at a subvector (u > of internal D.O.F. Then the total 

a 

displacement vector is partitioned into 


r > 


* u 

u 


c 



u 

1 


a 

u i J 


U 


where {u^} ™ subvector of internal D.O.F. where external loads 
are applied. 

{u } = complement of {u } in {u }. 

W a i 

In order to obtain good results for the set {u } where the 

a 

external loads are applied, the columns of the flexibility matrix 
corresponding to {u^}, in. the CMS model, must be complete. In 
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other words, the flexibility of the D.O.F. in {u > must be the 

a 

same in the CMS and FEM models. Similar* to free- interface 
components, the set of constraint modes [$ ] and fixed-interface 

c 

boundary flexibility vectors [Q^] do not supply the full 
flexibility for the D.O.F. in {u >. Thus residual attachment 

a 

modes have to be added to obtain the desired full flexibility. 
Hence, the contribution of the constraint modes [4 ] and the 

C 

fixed-interface boundary flexibility vectors [Q^] to the 

flexibility matrix, must be obtained. 

Let [T] be the matrix containing [$ ] and [Q ] in its columns 

c £ 


\ l *c ] 

[ QJ 1 = 

[I ] 

CC 

I 0 J*' 

l ° 

1 J 

[$,] 
L ic 

[Q]£ - 


The size of [T] is (n) x (c + l) . Where n is the total number of 
physical D.O.F. in component, c is the number of interface D.O.F. 
and l is equal to the number of the fixed- interface boundary 
flexibility vectors. Similar to equation (3.27), the contribution 
of [T] to the flexibility matrix is given by: 

-l 

tg k l = [T] [ [T] T tk] [T] j [T] T (5.5) 

where [g^l = contribution of [T] to the flexibility matrix. 
Substitute equation (5.4) in equation (5.5) and note that 

t®. 1 = ~ tk 1 1 [k ] and [k ] T = [k J 
Ic il ic ic cl 



the fixed- interface boundary flexibility vectors 




I 
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Let [g] be the total flexibility matrix of the component. In 
Appendix B, [g] is derived in terms of the stiffness submatrices 


of the fixed- interface component. 

Thus, 


£g i 


[gi = 

cc 



m . 


(5.9) 


where [g ] = [k ]-[k ][k J _1 [k ] 

cc L cc cl 11 ic 


-1 


[ <y ■ - [ny-tfcjnc.rikif tkjtk] 


-1 


[g ,i - tk ,r‘ * 


-i 


Hence the unrepresented flexibility matrix, or the residual 
flexibility matrix [g ] is given by: 

d 


[g a ] = tgj - [g„] 


r (o i 

CC 


to,] 

1 c 




[Q,][Q,] T 


(5. 10) 


For fixed-interface components, the set {u > is considered fixed. 

C 

Hence the residual flexibility matrix for the fixed- interface 
component will be [g^]^ 

t *a 1 „ * [k ,r - tQ t ItQ / ts.iu 

which is the full flexibility matrix for the fixed component, 
minus the contribution of the fixed- interface boundary flexibility 
vectors to the flexibility matrix* The residual flexibility 
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matrix given by equation (5.11), is used to obtain the residual 
attachment modes [# ] for the fixed component. This is achieved 
by applying unit forces at the subvector (u > of the physical 
0* F"" Hence the residual attachment modes [’p ] are the columns 

d 

of the residual flexibility matrix [g ] , corresponding to the 

D.O.F. in {u^}. It should be noted that the final displaced 

shapes formed by [* ], [QJ and [* J, will provide the full 
flexibility representation for the D.O.F. in {u } where the 

a 

external loads are applied. Thus the prediction of responses of 
D.O.F. in {u^} by using this formulation, is expected to be better 
than a formulation which does not use the residual attachment 
modes. If {p^} and {p^l the generalized coordinates 

associated with [Q^] and the sizes of {p^> and {p d > are it) 

and (a) respectively, where (£) is the number of boundary 
flexibility vectors in [Q^J, and (a) is the number of D.O.F. in 

{ u^} where the external loads are applied. The transformation of 

coordinates for the fixed- interface component, for case of 
external applied loads is given by: 


N 

► - m - 

u 

c 

p i 

► — 

ti ] 

cc 


[0 ] 1 
ca 


» « 
U 

c 

P# 

(5.12) 

L u i J 


. p d - 


[*,] 
L lc 

"V 

[* ] 
d J 


4 

l p J 



where [T] = transformation matrix of order (n) x (c + l + a). 
The transformation matrix [T], is used to transform the property 
matrices and the load subvector {f^}, as in equation (5.2), where 
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(5. 13) 


where {f^> = subvector of external loads applied at internal 

degrees of freedom {u }. 

a 

5 .3) A Comment on Wilsons Load-Dependent Ritz Vectors : 

As mentioned in chapter two, Wilson et al. [15] applied the 
load-dependent Ritz vectors to fixed-interface CMS. Only 
constraint modes were used for components having no external loads 
applied at internal D.O.F. For components having external loads 
<f > applied at a subset {u > of the internal D.O.F. {u }, 

a a i 

load-dependent Ritz vectors were generated from those applied 
loads. This section proves that the load-dependent Ritz vectors, 
the subsequent vectors generated from them and the constraint 
modes, will supply flexibility properties for the D.O.F. in {u >, 

a 

greater than that of the corresponding finite-element model. Thus 
higher displacement responses are obtained for the D.O.F. in (u }, 

a 

by using this model. 

Consider the same notation used in the previous section. 
Starting only with constraint modes [$ ], as displacement shapes 
for a component, we can proceed in a similar manner to that of the 
previous section, to obtain the contribution of the constraint 
modes to the total flexibility matrix. It was found to be equal 
to (compare to equation (5.8)): 
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[gj 

[g. ] , 

k cc 

k cl 

‘SlAc 

‘v,, 


where [ g ] 

lc cc 




(5. 14) 


‘S'- 


* ‘«.>L • 




-1 




= tk 




-i 


If [g] is the total flexibility given by equation (5.9), then the 
residual flexibility matrix [g^] (or* the unrepresented 
flexibility), is given by 


tg d ) = [gl - (g k l 


[0 ] 

CC 

to,] 

lc 



(5. 15) 


Which indicates that the constraint inodes supply full 

flexibility to the submatrices [g ] , [g ] and [g J, and an 

cc lc Cl 

incomplete flexibility to the submatrix [gj. Wilson et al. [15] 
used this matrix to generate the starting load-dependent 

vector from the external applied load. Thus the starting vectors 
will be the columns of the residual flexibility matrix, [g ] » 

d 

[kj \ corresponding to the D.O.F. in {u }, where the external 
applied loads act. Any subsequent vectors are generated according 
to Wilson’s algorithm given in chapter two. It should be noted 
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that the starting vectors, by themselves, will supply the full 
flexibility representation for the D.O.F. in {u > and ary 
subsequent generated vectors will add to the flexibility of the 
D.O.F. in {u^}. Thus the resulting component will have an "excess 
flexibility" for the D.O.F. in {u 1 and hence a higher response 

ft 

prediction for the D.O.F. where the external loads are applied. 

In order to obtain a component having a more realistic 
representation of the flexibility, another approach is suggested. 
This approach is similar to that taken in chapter three for 
free-interface components, in the boundary flexibility method of 
CMS. Assume, for simplicity, that the external load is applied 
only at D.O.F. number j from the set of degrees of freedom in 
{u^}. The static deflected shape obtained by this load is by 

definition equal to { , which is the attachment mode shape for 

D.O.F. number j. The { <p^} vector is equal to the column of the 

flexibility matrix, [g^] = [kj 1 , corresponding to D.O.F. number 
J* Is the starting load-dependent Ritz vector used by Wilson 

et al. [15]. In this formulation the starting vector will instead 
be 


{x*> = [kj" 1 [mj (5.16) 

where {x^ = starting load-dependent Ritz vector. 

It can be seen that {x*> is the inertial loading of the static 
response to the applied external load, and not the response 
itself, as in Wilson’s method. This is exactly similar to that 
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presented for free-interface components, in the boundary 
flexibility method of CIS (see section 3.3. 1.1.). The subsequent 
vectors are generated from the {x*> vector by using Wilson’s 
algorithm. It should be noted that orthogonal izat ions and 
normalizations are carried out with respect to the stiffness 
matrix [k^, instead of the mass matrix [mj. Thus a new simple 

residual flexibility matrix is obtained (as in the previous 
section) : 

[g„ 1 « tgj - [X][X] T 

d N d 

- [kj" 1 - [X][X] T (5.17) 

where tg^ ] = new residual flexibility matrix. 

[X] = matrix containing the generated load-dependent 

Ritz vectors in its columns. 

This new residual flexibility matrix is used to generate residual 
attachment modes [* J, which will be, in this case, the columns of 
the tg d J matrix corresponding to D.O.F. number J. The final 
transformation of coordinates, in its general form, is given by: 




f * 

U 


u 


c 


c 

• - [TJ - 

p * 

► = 

u 



i J 


p d J 



(I 1 [0], [0 ] 

cc cl ca 

[•] tx] [• ] 

1c d 


(5. 18) 


where (p^> = generalized coordinates associated with (£) 

generated load-dependent vectors in [X]. 

<P d } * generalized coordinates associated with (a) 

residual attachment modes [# ], corresponding to 

d 
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D.O. F. in {u^} where external loads are applied. 

[T] = transformation matrix of order (n) x (c + l + a). 

It should be noted that by this approach, we can ensure a full 
flexibility representation for the D.O.F. in {u >. The reason for 

a 

this is that the contribution of the generated load-dependent Ritz 
vectors to the flexibility matrix, was taken into consideration in 
equation (5.17). The property matrices and the external load 
vector are transformed, using the [T] matrix as before. 

5 .4) Sample Problem : 

The structure used to illustrate the forced vibration 
response by the boundary flexibility vector method of CMS, is the 
one shown in Figure (4.11), in the previous section. The external 
loading function applied is a step function, shown in Figure 
(5.1). The method used for integrating the equations of motions, 
is the Newmark step by step integration scheme (see Bathe [18]), 
with constants a and 5 equal to 0.25 and 0.5 respectively. This 
method is chosen because it is unconditionally stable for the 
given values of a and 5. The time step chosen for integration is 
0.02 seconds, which is less than the first period of the structure 
divided by ten. Proportional damping is used (see [16] and [18]) 
according to the following equation: 

[C] = a [M] + p [K] (5. 19) 

Where [C], [ M] and [K] are damping, mass and stiffness matrices 

respectively, for the whole system. a and £ are constants to be 
determined from two given damping ratios corresponding to two 
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unequal natural frequencies (see Appendix C). Two levels of 
damping are used; a low level corresponding to £ = 0.002 and 
^8 = 15 3X1(1 3 higher level corresponding to £ = 0.05 and 

? 8 = 0 . 15 . and £ are the damping ratios for natural 

frequencies one and eight respectively. For the above two levels 
of damping, the corresponding proportional constants are (a^ = 
— 7 • 489 , 0 ^ = 0 . 0124 ) and (a^= ” 3.584 , 0^ = 3-3033) respectively 
(see Appendix C). 

Three loading cases are considered, for each damping level. 
They are as follows: 

1 ) Load acting at node 31 in the global Z direction. Node 31 is 
an interface D.O.F. in components one and two. This case is 
solved by three methods; method one is by direct integration of 
the equations of motion, method two is by conventional CMS by 
using constraint modes and three fixed— Interface normal modes per 
component and method three is by the boundary flexibility method 
of CMS by using constraint modes and three fixed- interface 
boundary flexibility Ritz vectors. For all methods, Z-direction 
displacement responses are obtained at nodes 13 , 31 and 112 (see 
Figures ( 5 . 2 ), ( 5 . 3 ) and ( 5 . 4 ) for case of low damping, and 

Figures ( 5 . 5 ), ( 5 . 6 ) and ( 5 . 7 ) for case of high damping). The 

results obtained are almost the same for the three methods. For 
all responses, the case of high damping converges to the static 
solution after few cycles, while the case of low damping keeps 
oscillating about the static response, and requires a larger 
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number of response cycles in order to converge to It. For the two 
damping levels, the maximum and static responses, for each given 
node, obtained by the three different methods are shown in Tables 
(5. 1) and (5.2). 

2) Load acting at node 10 in the global Z-direction. This case is 
also solved by the above three methods. The Z-direction 
displacement responses are calculated at nodes 10, 84 and 112 (see 
Figures (5.8), (5.9) and (5.10) for low damping case and Figures 
(5.11), (5.12) and (5.13) for high damping case). Also the 
results obtained are almost the same for all three methods. The 
maximum and static responses obtained by each method, for both the 
two damping levels, are presented in Tables (5.3) and (5.4). 

3) Load acting at node 16 in the global Y-direction. This case is 
solved by the above three methods and by smother fourth method. 
The fourth method is the boundary flexibility method of CMS, with 
two boundary flexibility vectors, one residual attachment mode and 
constraint modes for component one where the load is applied, and 
three boundary flexibility vectors and constraint modes for each 
of the other three components. The Y-direction displacement 
responses are obtained at nodes 16, 87 and 99 (see Figures (5.14), 
(5.15) and (5.16) for low damping case and Figures (5.17), (5.18) 
and (5.19) for high dsunping case). Tables (5.5) and (5.6) show 
the maximum suid static responses obtained by all methods, for the 
two damping levels. For nodes 87 and 99, the results obtained by 
the four methods are almost the same, for the two levels of 
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damping. For the response at node 16, the results are not the 
same, as can be seen from Figure (5.14) and Table (5.5) for case 
of low damping, and Figure (5.17) and Table (5.6) for case of high 
damping. The results obtained, for the two levels of damping, by 
the conventional CMS method, indicates the maximum response is 
about 11% lower and the static response is about 14.5% lower than 
the direct integration results. Better results are obtained by 
the boundary flexibility method of CMS, in which three boundary 
flexibility vectors are used for each component. The results 
obtained by the boundary flexibility method of CMS, by using one 
residual attachment mode for component number one, are in the best 
agreement with the results obtained by the direct integration 


method. 
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In The 2-Direction (Low Damping) 



o.a 
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in The 2-Direction (Low Damping) 
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in The 2-Direction (Low Damping) 



Disp. 



Figure ( 5 . 5 ) - Response at Node 13 in The Z-Dlrection 
for The Case of Load Acting at Node 31 
In The Z-Directlon (High Damping) 



Disp. 



Figure (5.6) - Response at Node 31 in The Z-Direction 
for The Case of Load Acting at Node 31 
in The Z-Directlon (High Damping) 
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in The 2-Direction (High Damping) 




€00 


97 



•d»IQ 


in The 2-Direction (Low Damping) 
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for The Case of Load Acting at Node 10 
in The 2-Direction (Low Damping) 
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in The 2-Direction (Low Damping) 
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in The Z-Direction (High Damping) 
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for The Case of Load Acting at Node 10 
in The Z-Direction (High Damping) 



Conv. CMS 



time 


Figure (5.14) - Response at Node 16 in The Y-Direction 
for The Case of Load Acting at Node 16 
in The Y-Direction (Low Damping) 
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Conv. CMS 
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Figure (5.17) - Response at Node 16 in The Y-Direction 
for The Case of Load Acting at Node 16 
in The Y-Direction (High Damping) 
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Responses 

Method 

Direct 

Integration 

Conven. 

CMS 

B.F. CMS 

max. 

0. 19363 

0. 19382 

0. 19355 

static 

0. 10763 

0. 10763 

0. 10791 

Z 31 max. 

0.45955 

0.45959 

0.45967 

Z 31 static 

0.27631 

0.27631 

0.27631 

Z 112 max - 

0.06321 

0.06320 

0.06320 

Z 112 static 

0.02718 

0.02718 

0.02718 


Table (5.1 ) 


Maximum And Static Responses for 
Load at Node 31 in Z-Direction 
With Low Damping. 
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Method 

Responses 

Direct 

Integration 

Conven. 

CMS 

B.F. CMS 

Z 13 maX ' 

0. 18485 

0. 18503 

0. 18478 

Z ^ static 

0. 10699 

0. 10699 

0. 10699 

Z^ max. 

0.44217 

0.44223 

0.44231 

Z^ static 

0.27538 

0.27538 

0.27538 

Z 112 maX ' 

0.05895 

0.05894 

0.05894 

Z 112 static 

0.02670 

0.02670 

0.02670 


Table (5.2 ) - Maximum And Static Responses for 

Load at Node 31 in Z-Direction 
With High Damping 


















Responses 

Method 

Direct 

Integration 

Conven. 

CMS 

B.F. CMS 

max. 

0.08345 

0.08313 

0.08313 

static 

0.05265 

0.05236 

0.05236 

Z 84 max. 

0.06310 

0.06310 

0.06309 

Z g4 static 

0. 02722 

0.02721 

0.02721 

Z 112 maX ' 

0.02667 

0.02663 

0.02663 

Z^gStatic 

0.00192 

0.00192 

0.00192 


Table (5.3 ) — Maximum And Static Responses for 

Load at Node 10 in Z-Direction 
With Low Damping. 
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Method 

Responses 

Direct 

Integration 

Conven. 

CMS 

B.F. CMS 

Z 10 max. 

0.07661 

0.07625 

0.07624 

static 

0.05276 

0.05247 

0.05246 

Z 84 max - 

0.05905 

0.05905 

0.05905 

Z 84 static 

0.02715 

0.02715 

0.02715 

Z 112 max - 

0.01931 

0.01928 

0.01928 

Z 112 static 

0.00184 

0.00184 

0.00184 


Table (5.4 ) - Maximum And Static Responses for 

Load at Node 10 in Z-Direction 
With High Damping. 
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Responses 

Method 

Direct 

Integ. 

Conven. 

CMS 

B.F. CMS 

B.F. CMS With 
Resid. Mode 

i 

B 

CO 

rH 

>• 

0. 08908 

0.07925 

0.08695 

0.08781 

Y^g static 

0.06522 

0.05583 

0.06287 

0.06386 

Y g7 max. 

-0.02239 

-0.02239 

-0.02251 

-0.02247 

Yg 7 static 

-0.00353 

-0.00353 

-0.00353 

-0.00353 

Y 99 max - 

-0.01219 

-0.01212 

-0.01231 

-0.01229 

Ygg static 

-0.00193 

-0.00193 

-0.00193 

-0.00193 


Table (5.5 ) 


Maximum And Static Responses for 
Load at Node 16 in Y-Direction 
With Low Damping. 













Method 

Responses 

Direct 

Integ. 

Conven. 

CMS 

B.F. CMS 

B.F. CMS With 
Res id. Mode 

Y„ max. 
lb 

0.08793 

0.07812 

0.08578 

0. 08665 

static 

0.06555 

0.05611 

_ 

0.06318 

0.06417 

Y 87 max. 

-0.02161 

-0.02159 

-0.02172 

-0.02168 

Yq 7 static 

-0.00351 

-0.00351 

-0.00351 

-0.00351 

Y 99 maX ’ 

-0.01048 

-0.01043 

-0.01060 

-0.01058 

Ygg static 

-0.00194 

-0.00194 

-0.00194 

-0.00194 


Table (5.6 ) - Maximum And Static Responses For 

Load at Node 16 in Y-Direction 
With High Damping. 























Chapter (6) 


S ummary and Conclusio ns 

A new method of component mode synthesis was presented in 
this work. The new method, namely the boundary flexibility vector 
of CMS, is based upon a set of static Ritz vectors (boundary 
flexibility vectors) as generalized displacement shapes for 
components. The generation of these vectors does not require the 
solution of the eigenproblem associated with the component, as in 
the case of conventional methods of CMS. The formulation of the 
new method for the free vibration problem as well as the forced 
vibration problem, was presented for both fixed and free-interface 
components. A comparison of the number of operations required to 
obtain the boundary flexibility vectors versus normal modes of 
vibration, was presented. The comparison showed a substantial 
reduction of the number of required operations by using the 
boundary flexibility vectors, instead of eigenvectors in CMS. 

Numerical examples for the free vibration problem were 
presented in chapter four. The examples were solved by the 
boundary flexibility method of CMS and by other different methods. 
Results indicated that by using the boundary flexibility method, 
more accurate results could be obtained with a substantial saving 
in CPU time. The new method was applied to a substantially large 
structure, in chapter four. The results obtained by the boundary 
flexibility method of CMS, were in good agreement with 
finite-element and conventional methods of CMS. The saving in CPU 
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time gained by using the boundary flexibility vectors, instead of 
eigenvectors, was computed for this problem. An average saving 
of CPU time between 75% and 80% could be attained, for generating 
the same number of vectors. 

A numerical example was presented for the forced vibration 
problem. The results obtained by the boundary flexibility method 
of CHS, were in good agreement with the results obtained from the 
direct integration of the equations of motion. In all cases 
considered, the results obtained by the boundary flexibility 
method of CMS were almost the same and sometimes superior to those 
obtained by conventional methods of CMS. 

It was proved by this work, that the boundary flexibility 
method of CMS could be applied to free and forced vibration 
problems. It was shown that accurate results could be obtained by 
applying the method, with a substantial saving in CPU time 
compared to conventional methods of CMS. 

A general outline for continuation of this work may be the 
following: 

1- Develop a criteria to specify which of the constraint or 
attachment modes, are to be chosen to generate the boundary 
flexibility vectors. This criteria should give an indication of 
how much a constraint or an attachment mode will contribute to the 
response. 

2 Investigate methods for refining the generated boundary 
flexibility vectors. One idea suggested is to investigate the 
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generation of the vectors from actual interface stiffness 
properties between components, instead of assuming the interface 
either fixed or free. 
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Appendix A 


Comparison of The Number of Operations For Generating 
The Boundary Flexibility Vectors Versus Eigenvectors 

The number of operations are derived using some basic 
criteria in reference [19]. 

A . 1) Extraction of k Eigenpairs : 

For the deteminant search and inverse iteration methods, see 
Bathe [18] 

Assume: n = total number of degrees of freedom 
k = number of required eigenpairs 

y = number of iterations in determinant search method 
required to obtain eigenvalue i 
= number of inverse iterations required to obtain 
eigenvector i 
[K]= stiffness matrix 
[M]= mass matrix 
{Y^}^ eigenvector number i 

p = amount of shift applied to stiffness matrix 
The following steps are required to obtain each eigenpair, thus 


they 

Step 

are repeated k times. 
Operat ion 

No. of Mult, or 
Divisions 

No. of Add. or 
Subtractions 

1 

Apply shift 




[Kj = [ K ] - p [M] 

2 

n 

2 

n 


Triangularize [KJ 

2n 3 + 3n 2 - 5n 
6 

3 

n - n 
3 


Total for step 1 

2n 3 + 9n 2 - 5n 

n 3 + 3n 2 - n 



6 

3 
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Step 


Operat Ion 


Step 1 is repeated from 
j = l -...y. 


Obtai 


n [Kj 


l 

-l 


Compute {X> +i = [KJ -1 {Y } 


Compute {Y> = [M] {X} 

j+i J+i 


Compute p 


{X iI: {Y 5 


J + 1 {X> T {Y> 

J+i j+l 


m 


Compute {Y} = 
J + i 


j + i 


1/2 


,x i < y s .. 


Total for step 3 

Step 3 is repeated from 
j = 1 x 


No. of Mult, or! 
Divisions 


. 3 

4n - n 


n 


n 


2n + 1 


2n 


2n + 4n + 1 


No. of Add. or 
Subtractions 


8n 3 - 9n 2 + n 


2 

n - n 


2 

n - n 


2n - 2 


n - 1 


2n + n - 3 


Then, total number of multiplications and divisions = 

. s. r y. 

T 


k r y -i 

(4n -n) + £ (2n 3 +9n 2 -5n) + x i (2n 2 +4n +1) j 


(A. 1) 


Total number of additions and subtractions = 

* r y. 


(8n 3 -9n 2 +n) + £ (n 3 +3n 2 -n) + x t (2n 2 +n- 3) j 


(A. 2) 
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A . 2) Generation of l Boundary Flexibility Vectors : 
Assume { > = boundary flexibility vector number i 


Step Operation 

1 Obtain [K]* 1 

2 Compute {X^ = [M] {q^} 

Compute {q^ = (K] 

Normal izing 

• * T 

GM = {q*>* [M] {q*} 

i + 1 i + 1 

<V.,- Cm { <U 

Total for step 2 

Step 2 is repeated for 
i = 1 I 

3 Orthogonal izing 

Compute c^= {q^> [M]{q^> 

• » • 

Compute {q } 

Total for step 3 

Where j = 1 i-1 

Step 3 is repeated for 
i = 2 I 


No. of Mult, or 

No. of Add. or 

Divisions 

Subtractions 

. 3 

4n - n 

8n 3 - 9n 2 + n 

3 

6 

2 

2 

n 

n - n 

2 

2 

n 

n - n 

2 

2 , 

n + n 

n - 1 

n 


3n 2 + 2n 

3n 2 - 2n - 1 

2 

2 . 

n + n 

n - 1 

n 

n 

2 ~ 
n + 2n 

2 

n + n - 1 
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Total number of repetitions for step number 3 = (1+2+3. 


£-i 

■ I k 


k= 1 


l U- 1 ) 


Hence total number of multiplications and divisions = 

t ) n 2 ♦ ( gil- V* - , „ 


Total number of additions and subtractions = 


( 4- ) n 3 + ( il^-1 ) n 2 + ( 3^151^1 , n _ { , 


U-U) 


(A. 3) 


2 


6 


2 


(A. 4) 




Substitute in equation (B.4), then 





o 
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Appendix C 


Proportional Damping Matrix 

Proportional damping matrix [C], (see Bathe [18] and Clough 
and Penzien [16]), is given by 

[C] = a [M] + /3 [K] (C. 1) 

where [M] = mass matrix 

[K] = stiffness matrix 

a and /3 are constants to be determined from two given 
damping ratios corresponding to two unequal frequencies of 
vibrations. 

Assume mode shape {<p^t has a frequency Pre and post multiply 

both sides of equation (C. 1) by and respectively 

2 u> l $ l = (a [M] +0 [K] ) 

2 = a + £ cf (C. 2) 

where (; = damping ratio of mode shape i. 

By substituting, in equation (C.2), any two unequal frequencies 
and their corresponding assumed damping ratios, a and £ can be 
calculated from the resulting two simultaneous equations. The 
damping ratio of any other natural frequency is obtained from: 

a + /3 ui 2 

= - (C. 3) 

2 u 

i 
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